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Abstract 

We develop a dynamical theory, based on a system of ordinary differential equations describing the 
motion of particles which reproduces the results of quantum mechanics. The system generalizes the Hamil- 
ton equations of classical mechanics to the quantum domain, and turns into them in the classical limit 
h —¥ 0. The particles' motions are completely determined by the initial conditions. In this theory, the wave 
function ip of quantum mechanics is equal to the exponent of an action function, obtained by integrating 
some Lagrangian function along particle trajectories, described by equations of motion. Consequently, the 
equation for the logarithm of a wave function is related to the equations of motion in the same way as 
the Hamilton- Jacobi equation is related to the Hamilton equations in classical mechanics. We demonstrate 
that the probability density of particles, moving according to these equations, should be given by a stan- 
dard quantum-mechanical relation, p — l?^ 2 . The theory of quantum measurements is presented, and the 
mechanism of nonlocal correlations between results of distant measurements with entangled particles is 
revealed. In the last section, we extend the theory to particles with nonzero spin. 
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1. Introduction 



According to standard quantum mechanics!!! (QM), the state of every physical system is described by a wave 
function, whose time evolution is determined by the Schrodinger equation. In this paper we will consider only 
closed systems, for which the description by a wave function is sufficient. We know how to set up experiments 
with a known initial wave function, and then using the Schrodinger equation we can calculate it at any later 
moment. However, contrary to, say, an electric held in an electromagnetic wave, or a held of pressure in a sound 
wave, the wave function is not an object of observation and measurement. Consequently, besides Schrodinger's 
equation, the theory additionally includes a set of rules, specifying the results of experiments with quantum 
systems in terms of their wave functions. These rules were developed in late 1920-s and collectively named the 
(statistical) interpretation of QM. Thus, the theory has two parts: Schrodinger's equation and interpretation. 

Such structure of the theory may be viewed in various ways. The standard attitude consists of the faith 
that the described construction constitutes the desired complete and fundamental law of nature. However, 
there are a number of objections that may be raised against this point of view: 

- It seems natural to expect from a fundamental theory that it reflects all observable elements of physical 
reality and gives the law of evolution for them. Thus the very fact that QM is formulated in terms of wave 
functions, which cannot be directly observed, and requires an interpretation that establishes a connection 
between wave functions and results of experiments creates doubts in its fundamental character. 

- This interpretation is a separate and independent part of the theory's foundation, whereas it seems desirable 
for a fundamental theory to allow the derivation of all its experimental consequences by pure math from the 
dynamical laws of evolution alone. 

- By necessity the interpretation, which describes the response of an approximately classical apparatus to its 
interaction with a quantum system, is expressed in classical terms. However the behavior of any apparatus, 
which is just a physical object built up of atoms, should be derivable from QM. Consequently, QM contains an 
unacceptable for a fundamental theory logical vicious circle: in the words of a classic textbook pQ: "... quantum 
mechanics occupies a very unusual place among physical theories: it contains classical mechanics as a limiting 
case, yet at the same time it requires this limiting case for its own formulation." 

- The interpretation happens to be probabilistic, thus employing a series of similar experiments, possibly per- 
formed in different places and at different times, to establish the meaning of a wave function in the experiment 
at hand. The wave function in this particular experiment, however, certainly appears relevant. It therefore 
seems desirable for the theory to define the meaning of a wave function in every individual experiment, without 
reference to its repetitions (especially when such repetitions are clearly impossible, such as when discussing the 
wave function of the universe) which QM fails to do. The fact that a wave function must have a nonstatistical 
interpretation in "internal" terms also clearly follows from the utility of a concept of a wave function of quarks 
conhned inside hadrons. 

- According to this interpretation, during a measurement the wave function abandons the unitary law of evo- 
lution, which it normally follows, and suffers a collapse. However, the conditions under which this change of a 
character of evolution happens, are not specified in QM, and attempts to formulate such conditions have not 
been convincingly successful. 

- The collapse occurs randomly into different possible states, but QM does not explain the reason for this 
randomness. Consequently, the values of corresponding probabilities (which one would expect to see among 

1 In this article we will use, in order of their appearance, the following abbreviations: QM — quantum mechanics, ODE — 
ordinary differential equation, PDE — partial differential equation, QHJE — quantum Hamilton-Jacobi equation, "PAQD — 
analytical quantum dynamics in infinite phase space, OSFI — one-step Feynman integral, CD — Cartan distribution, HC — 
Hamiltonian conditions, DBBT — de Broglie-Bohm theory, FDS — full description space, RDS — reduced description space. 
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the results of a fundamental theory) in QM are not derived, but postulated, or, in other words, are taken from 
experiment. 

- According to QM, unless a system is in an eigenstatc of a measured quantity, the result of its measurement 
does not exist before the measurement is done, but is rather created during the measurement. In some cases 
this is completely obvious — see a discussion of spin measurement for a particle with spin 3/2 in [2J. This 
means that a measurement is not a fundamental unanalyzable primitive, but a nontrivial physical process for 
which QM fails to give an adequate description. 

- Moreover, this process of measurement produces nonlocally correlated results for measurements performed 
with space-like separated entangled particles, but QM does not describe any mechanism which causes these 
correlations. 

Thus it appears that the standard combination of Schrodinger's equation and statistical interpretation is 
too complicated, artificial and, in the words of John Bell "unprofessionally vague and ambiguous" [3J, while the 
nature obviously prefers simple, natural, and clear fundamental laws. Consequently, we suggest in this paper to 
regard the situation in the following alternative way: The statistical interpretation does not, of course, follow 
from the Schrodinger equation (simply because the latter only describes the behavior of a wave function), 
but is a generalization of results of observations and experiments. Our trust in statistical interpretation is 
based on its agreement with experiment, and only on this agreement. Therefore, the interpretation is a 
phenomenological part of quantum theory, and so the whole existing theory is semi-phenomenological. Then 
to this semi-phenomenological theory the above objections are inapplicable, while at the same time there 
remains a possibility that the nature is ruled by the other, "simple, natural, and clear" fundamental theory, 
from which statistical interpretation (and maybe Schrodinger's equation as well) follows. 

The conclusion about the phenomenological (or, as it is often called, "pragmatic" [4J nature of exist- 
ing QM may also be drawn from the works devoted to its foundations. We read, for example, in Bohr [5]: 
"Strictly speaking, the mathematical formalism of quantum mechanics and electrodynamics merely offers rules 
of calculation for the deduction of expectations about observations obtained under well-defined experimental 
conditions specified by classical physical concepts", or, in a frequently quoted more recent paper [BJ: "...quan- 
tum theory does not describe physical reality. What it does is provide an algorithm for computing probabilities 
for the macroscopic events", "...the time dependence of the wavefunction does not represent the evolution of a 
physical system. It only gives the evolution of our probabilities for the outcomes of potential experiments on 
that system." Thus according to these works, QM describes the results of our observations of electrons, atoms, 
etc. Which theory, then, describes these particles, which we observe, themselves! 

Of course, neither these quotations, nor the arguments presented above, can prove that a better theory 
is needed. They can, however, motivate a search for such a theory. Indeed, it is hard to help feeling that 
peculiar features of QM are the consequences of a fact that it misses some important part of a complete 
theory, a part which is substituted by a phenomenological description of the way it works. In this searched-for 
complete theory, the wave function must have a definite meaning in every individual experiment, and the 
theory must explain the nature of randomness and derive the standard quantum-mechanical expression for 
probability. The expression for probability will thus become just a property of a wave function, rather than a 
basis for its interpretation. This theory should also explain the properties of quantum measurements, describe 
the mechanism which creates nonlocal correlations, and fill all other gaps listed above; in particular, it should 
contain an image of every observable element of reality and predict its behavior directly from the theory's 
dynamical laws, without the need for any special interpretation. Compared to QM, such a theory would be 
much less vulnerable to suspicions of being a mere semi-phenomenology, and it is a goal of this paper to present 
a theory which appears to satisfy these demands. 
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Before discussing this theory, we recall the Hamilton- Jacobi equation 

% + (VP) 2 + U = (1.1) 
at 2m 

for an action function p(x, t) in classical mechanics. This function does not describe any individual trajectory 
and motion of a particle along it; rather, it describes a family of such trajectories, of which none can be singled 
out given an action function alone. Individual trajectories and particles' motions along them are described in 
classical mechanics by Hamilton ordinary differential equations (ODEs). Given these trajectories, the action 
function, which solves the Hamilton- Jacobi equation, may be obtained by integrating a Lagrangian function 
along them. On the other hand, the trajectories of a family, described by an action function p(x, t), may be 
reconstructed from this function using the equality p(x, t) — Vj>(x, t), which says that the momentum p at 
point (x, t) of any trajectory is equal to the gradient of an action function at this point. Thus the Hamilton 
and Hamilton- Jacobi equations represent two parts of the same theory — classical mechanics — the former 
describing particle motion along classical trajectories, and the latter, properties of the families or ensembles 
of these trajectories. 

Now we formulate the basic idea of the present approach. Its initial step is purely mathematical. Namely, 
it is shown in a second section of the paper, that similar to the case of first order partial differential equations 
(PDEs) with one unknown function, such as just discussed Hamilton- Jacobi equation, for a large class of 
PDEs of second and higher orders the solution of equation may be represented as an action function, i.e. the 
value of the function p(x, t) that solves the equation may in every point be obtained as an integral from some 
"Lagrangian" function along the curve that leads to this point and is completely and uniquely determined by 
some system of ODEs. As is well known [H [9l [TD] , for PDEs of higher than first order the system of ODEs 
with such properties neither exists in the usual phase space with coordinates t,x l ,pi, where pi = dp/dx l are 
first derivatives of an unknown function, nor even in the same space extended by adding to its coordinates 
the derivatives of an unknown function up to any finite order. Such a system, however, exists in an infinite 
phase space, the coordinates of which include all possible partial derivatives of an unknown function, and a 
corresponding mathematical theory is developed in section 2. Although the very possibility of solving higher- 
order PDEs in this way was known for quite some time QT] , the specific form of solution presented in section 2 
seems to be new. In spite of the presence of an infinite number of variables and equations in the theory, it 
happens to be quite transparent and manageable; in fact, the theory is remarkably similar to the Hamiltonian 
formalism in classical mechanics and reproduces all its essential features. The theory of first order PDEs also 
can be formulated in an infinite phase space and turns out to be a special case of our theory, but in this case 
ODEs for x % and pi, i.e. for coordinates in the usual phase space, decouple from other equations and can be 
considered independently. Thus, we obtain a general Hamiltonian formalism that covers a large class of PDEs 
of first as well as higher orders on equal grounds. 

Returning now to physics consider, along with a wave function ip, a function p — (h/i) lnip. Clearly, this 
function contains the same information as ip, and may be used instead of it in all discussions. For a spinless 
particle of mass m in external potential U we have from the Schrodinger equation the following PDE for p(x, t): 

t + ^-(Vp) 2 + C/+^-Ap = 0. (1.2) 
at 2m 2im 

Except for the last term, proportional to h, this is the Hamilton- Jacobi equation for an action function in 
classical mechanics. On the other hand, this equation happens to belong to the class of PDEs considered 
in section 2, which have a solution in the form of an action function. In this situation the following main 
idea of the present approach emerges with an absolute inevitability: consider Eq. (11.21) as an equation for the 
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Equations of motion 
(ODEs) 



Equation for action function 
(PDE) 



Classical theory 
Quantum theory 



Newton or Hamilton equations 



Hamilton-Jacobi equation 



"PAQD equations of motion 



Quantum Hamilton-Jacobi equation 



Figure 1: The structure of classical mechanics and PAQD. The relations between equations of motion and 
equations for action functions are the same for both theories, and in the limit h — ¥ quantum equations turn 
into classical ones in both columns. In standard QM the lower left rectangle is absent and substituted by the 
statistical interpretation. 

action function in a new, quantum, theory, the wave function — as an exponent of the new action function 
(multiplied by i/h), the curves along which the Lagrangian function should be integrated to produce the 
action function — as particle's trajectories in the new theory, and ODEs that determine these curves — as 
new equations of motion, which correct Hamilton's equations. In exact analogy with classical mechanics, the 
resulting theory will have two sides: ODE side, represented by the equations of motion of particles along their 
trajectories, and PDE one, represented by Eq. (11.21) for the action function that describes, along with the 
wave function, ensembles of trajectories. In the following, Eq. (|1.2[) will be called "quantum Hamilton-Jacobi 
equation" (QHJE). The motion of particles takes place in an infinite phase space, V, defined in section 2, and 
the theory will be called analytical quantum dynamics in infinite phase space, or 'PAQD. The general structure 
of 'PAQD and classical mechanics is presented in Fig. [1] 

In view of described above new mathematical possibility of dealing with Eq. (jl.2p . it seems difficult to 
dispute that this approach certainly appears quite natural; one could say that by the very form of Eq. (jl.2p 
the nature pushes us in this direction. Had the possibility of solving second order PDEs in the way described 
above been known in 1926-1927, it is hard to doubt that this work would be done right then! Further, being 
based on equations of motion and only on them, the theory should be considered simple; for the same reason 
it promises to be clear and unambiguous. Thus, it seems worth the efforts to investigate the possibilities 
which may open in this direction; in doing so we will also finish the job left unfinished eighty five years ago 
due to such historical accident as an absence of a proper mathematical formalism at that time. Last but not 
least, we note that the theory is completely fixed by Eq. (|1.2[) for its action function (that is — fixed by the 
Schrodinger equation) and doesn't contain any additional freedom to improve its agreement with the second 
part of standard QM, the statistical interpretation. Therefore the fact that, as we will soon see, such agreement 
is nevertheless achieved (or, in other words, that statistical interpretation is deduced from the present theory) 
should be considered as a weighty argument in the theory's favor. 

The equations of motion of the theory are explicitly written down in the beginning of section 3, and the rest 
of the paper is devoted to demonstrating that QM may be understood as a theory of particles moving according 
to these equations, the difference between classical and quantum mechanics being the result of the different 
form and number of equations in these theories. Equations of motion arc followed in section 3 by their general 
discussion. As in classical mechanics, these equations are self-sufficient: given initial conditions, they define 
the particle's motion unambiguously, without any need for using a wave function or the Schrodinger equation. 
In the classical limit h — > 0, the equations of motion turn directly into the Hamilton equations of classical 
mechanics. The equations and particle trajectories live in an infinite phase space. We show that projections of 
these trajectories to physical space coincide with "Bohmian trajectories," introduced by de Broglie and Bohm 
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[T2l [T3] on very different grounds, and discuss the reiation between the present theory and that of de Broglic 
and Bohm. 

The class of PDEs covered by the theory in section 2 is very large, and so the question arises: what singles 
Eq. (II. 2[) out of this class? It is shown in section 4, that equations of motion may be obtained from a "one-step 
Feynman integral" (OSFI), combined with an appropriate variational principle. OSFI provides, therefore, an 
alternative starting point of the theory, bypassing the geometric theory of section 2. Moreover, OSFI is a 
functional of a Lagrangian function, which is necessarily classical, i.e., depends on the position and velocity of 
a particle only. Consequently, OSFI may be considered as a general source of quantum theories, obtained by 
"quantization" of corresponding classical theories, represented by Lagrangian functions. As further discussed 
in section 4, the theories, obtained in this way, will automatically exhibit familiar features of QM: superposition 
principle, path-integral representation, and wave-particle duality. We note that in a mathematical derivation 
of the latter feature, an infinite number of variables and equations in our theory, which initially appears to be 
a theory's disadvantage, plays a crucial role. 

Another consequence of OSFI is that a corresponding PDE may be obtained from a variational principle. 
This is shown in the beginning of section 5.1. By Noether's theorem, the symmetries of such a PDE lead to 
conservation laws. We then use a fundamental invariance of all PDEs, considered in section 2, with respect to 
a shift of the unknown function by a constant to derive a continuity equation. In section 5.1 this is done for 
a standard Hamiltonian of the Schrodinger equation, and in section 5.2 — in a general case, without using an 
explicit form of a Hamiltonian. In section 5.3, we use the current conservation to prove that a form |^| 2 dy 
is an integral invariant of our equations of motion, which replaces the canonical integral invariant (Liouville 
measure) d?p dV of the Hamilton equations in classical mechanics. 

Using the invariance of the form l^pdU, section 6 demonstrates that a probability density in configuration 
space should be equal to |-0| 2 . We give two proofs, the second one using the maximization of a specially 
introduced functional of probability density, analogous to the Gibbs entropy. We compare the situation in 
QM to the one in classical statistics. A brief review of equilibrium and nonequilibrium classical statistics is 
presented in the Appendix in a form convenient for such comparison. It is shown there that the repetition 
of steps which led to the expression of Gibbs entropy in QM leads in classical statistics (where for invariant 
measure one uses the form d 3 pdV, rather than |^| 2 c?U) to its standard classical expression. The probability 
density \ip\ 2 that maximizes Gibbs entropy in QM has, therefore, a status identical to that of a microcanonical 
distribution, which maximizes Gibbs entropy in classical statistics. Note that the difference between these 
distributions results from the difference between corresponding invariant measures, which, in turn, follows 
from the difference in equations of motion. Regarding the claims [T2HH] that the \ip\ 2 distribution in QM may 
arise in a way similar to relaxation to statistical equilibrium in classical statistics, the Appendix also shows 
that this relaxation is related to the growth of Boltzmann, rather than Gibbs, entropy, and is caused by the 
properties of macroscopic systems which cannot have any analogs in a one-particle theory. 

In section 7, the one-particle theory of the previous sections is generalized to multiparticle systems. We also 
discuss how the standard physical picture of quantum particles in a potential created by classical macroscopic 
objects emerges from our theory. 

Section 8 considers the theory of quantum measurements. Von Neumann's measurements with discrete and 
continuous spectra are considered in sections 8.1 and 8.2. The theory discussed there is a "PAQD-adaptation 
of the theory developed by Bohm [T^l [T3]. In section 8.3, the measurement of a particle's position by a 
photographic plate or in a bubble chamber, which is not a von Neumann's measurement, is considered, and 
its properties are discussed. In the end of this section we analyze the double-slit experiment discussed by 
Feynman [T5| and compare its results with "PAQD predictions. 
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Section 9 considers the mechanism of nonlocal correlations between the results of measurements, performed 
with space-like separated, but entangled, particles. We argue that the relativistic version of PAQD, although 
nonlocal, will be Lorentz invariant. 

Section 10 considers particles with spin. We show that their theory, which adequately generalizes the 
theory of spinless particles, may be developed based on extended configuration space, which includes, besides 
the particle's space position, also its "internal" SU(2) coordinates. 

Finally in Conclusion, we give a brief review of our theory, compare it with standard QM, and finish with 
several general remarks. 

2. Hamiltonian flow in infinite jet space 
2.1. Basic definitions and notation 

In this section we discuss the question of when the solution of a PDE system may be obtained, as in the case 
of the Hamilton- Jacobi equation, via solving some related system of ODEs. Consider evolutionary PDEs of 
the form 

-j^ + H r = 0, r = l,...,m, (2.1) 

where the p r are m unknown functions ("dependent variables") of n space variables q , . . . , q n combined into a 
vector q, and time t ("independent variables"), and H r are functions of t, q, and partial derivatives of unknown 
functions with respect to space variables up to some finite order. Denote these derivatives by corresponding 
multi-indices, as in p r ^ = d 2 p r /dq z dq 3 , and include in the set of all possible multi-indices an empty one, denoted 
as 0, which will correspond to the function p r itself. Use k for space indices, running from 1 to n, use r and 
s for function indices, running from 1 to m, and use Greek letters for multi-indices. The order of indices in a 
multi-index is arbitrary, and two multi-indices which differ only by permutation are considered to be the same. 
Correspondingly, only one such multi-index will be assumed to be included in a summation over all possible 
multi-indices. If a = ...%}■, let ai or ia be the "extended" multi-index i\i2 ■ ■ ■ ikh and if fx = jx . . . ji, let 
(T/i or [ia be the multi-index i\ . . . . . . ji- Let Cj, i = 1, . . . , n, denote the number of times index i is found 

in the multi-index a, so that a may be represented as u\ ones, followed by oi twos, etc. It is useful to think 
of the multi- index a as an n-dimensional vector with nonnegative integer components <7j. Summation over 
all possible multi- indices a then reduces to summation over all af. = a =o ■ l°1 denote the 

total number of indices in multi-index a, so \a\ = Y^i=i a i- = — d/dt. For every multi-index er, let 

cr! = Yli=i °i- an d ®cr — n"=i(^/^9 l ) <Ti - For any n-dimensional vector x, let x CT = nlLi( a; *) CT< - 

Let {p} denote the set of all unknown functions and all their derivatives. By analogy with classical 
mechanics, functions p r will be called action functions, or just actions, and their derivatives, momentums. 
Denote the set of all momentums, i.e., all p r a with cr ^ 0, by p, so H r = H r (t, q,p). Denote the space of 
independent variables q and t ("base space") by M, and the space of vectors q alone ("configuration space") 
by Q. Let J°° ("infinite jet space") be the space with coordinates t, q, {p}, i.e., all independent as well as 
dependent variables and all their space derivatives. Call the similar space V with coordinates t, q, p "infinite 
phase space." We assume that H r depends analytically on its arguments, and consider only analytic or real- 
analytic solutions of Eq. (|2.1j) . A mathematically rigorous treatment of geometry of analytic jets in J°° may 
be found in [TB]. Denote by J k a jet space of fc-times continuously differentiable functions, which includes 
derivatives only up to fc-th order ( "k-jets space" ) [7J |H1 [H] ■ We will also use the notation when it is necessary 
to indicate the number m of unknown functions explicitly, and when m — 1 we will drop the function number 
superscript in equations. We will consider every solution of Eq. (|2.1[) , with all its space partial derivatives, as 



creating a graph in J 00 . Denote such a graph by T and note it is an n + 1 - dimensional surface in J°° . We 
want to find out when such graphs can be usefully considered as formed by a congruence of curves described 
by a system of ODEs. 

2.2. Curves in the graph of a PDE solution 

It is easy to write an equation for an arbitrary curve which lies in the graph. Let q(i) be the projection of the 
curve to the base space M. Consider the operator of total differentiation 

r— 1 a ° 

where the second summation runs over all possible multi- indices a. At every point of J°°, the operator Di 
raises the partial derivative djdq % to a graph of an analytic function of q, which passes through the point 
(see [7]). In other words, if O is such a graph, {pe(q)} is a set of values of unknown functions and their 
derivatives at a point of G with a base coordinate q, and F(t, q, {p}) is some function in J°°, then 

AF(i,q,{ Pe (q)}) = AF(t,q,M)| M={pe(q)} . (2.3) 

By consecutive differentiation of Eq. (|2.1[) with respect to space variables we now obtain equations ( "prolongations'' 
of (I2.ip ) which describe the behavior of space derivatives p r a of the solution 

^ + K = 0, (2.4) 

where H r a — D a H r and D a is repeated total differentiation, i.e., D a — Y\^ =1 {Di) ai . Eq. (|2.4|) gives the value of 
the partial derivative of p r a with respect to t, while its partial derivative with respect to q l is, by definition, p r ai . 
Consequently, the time dependence of the T-image of the point q(t) on the base is described by the system of 
ODEf§ 

V ° dq* q + dt (2 .5) 

= vlitf-K, 

where by dot we denote the total derivative of a corresponding value with respect to t along the curve q(t) 
(see Fig.EJ). 

Given the PDE (|2.1jl . we constructed the ODE system (|2.5|) . which defines the evolution of values of 
unknown functions and their derivatives at a point q(i) moving (in an arbitrary way) through the base. It is 
easy to show that there also exists an inverse correspondence: if there is a point q(t) moving through the base, 
and a sequence of functions of time p%(t), which satisfies the system of ODEs (|2.5[) . then functions p s (r,t), 
which are the sums of Taylor series in a point q(t) with coefficients p%if), will satisfy PDE (|2.1j) . Indeed, these 
Taylor series are 

p s ^ =E^w( r -qwr- ( 2 - 6 ) 

Taking the time derivative and using Eq. (|2.5p for , one obtains, after some cancellations, 

^{v,t) = -^^MWll^wjhqfC, (2-7) 



2 Summation over repeated indices and multi-indices is assumed here and below unless stated otherwise. The summation will 
not be assumed when the expression with repeated indices or multi-indices stands on one side of an equality if the other side of 
the equality also has the same (multi-)indices used only once. 
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Figure 2: A graph T of solution of PDE (|2.1[) in the infinite jet space. When a point moves in the configuration 
space, the set of its coordinates q and values {p} of solution and solutions's derivatives at q moves in T according 
to the system of ODEs (|2~5j) . 



which by Eq. (|2.3[) and analyticity of H s is equal to —H s (t,r, {pp(r)}), where T is the graph of functions (|2.6j) 
in J°°, as required. 

We see that there exists a simple and general relation between the evolution of analytic functions, described 
by PDE (I2.ip . on one hand, and the evolution of coordinates q(i) and momentums p r a {t), described by ODE 
system (I2.5p . on the other. This relation is just a straightforward consequence of the structure of Eqs. (|2.ip 
and (|2.5p . and is satisfied for an arbitrary curve q(t). However, if one adds to system (|2.5p an additional 
equation, expressing q % through other variables, then this extended system of ODEs will be closed with respect 
to the evolution of all variables involved, and may be viewed as representing equations of motion that are 
generated by "Hamiltonians" H r in the same way as Hamilton equations are generated by a Hamiltonian in 
classical mechanics. Consequently, to obtain potential equations of motion, some plausible condition capable 
of providing such an equation is needed. 



2.3. Generalized Hamiltonian fields and equations of motion 

To formulate the condition described above, we will use concepts and notation of the geometric theory of 
PDEsi [3 M UM EH OS- The set {A, d/dpl, d/dt} is a basis for the tangent bundle T( J°°). Define the 
sequence of 1-forms 

Z r a = dp r a - p r al dq\ (2.8) 
The set {dq\ u£, dt} is a basis for the cotangent bundle T*(J°°). The duality relations 
Di^dqi = Sij, £><_lD£ = 0, D.Jdt = 0, 

= 0, ^_i«3J = 0, |_.di = 1, 



3 A reader unfamiliar with this theory may skip directly to Eq. 1 12.2211 , which is a desired expression for velocity. An elementary 
justification of this expression is presented in a footnote immediately after it. 
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then follow at every point of J°°. In the following, let T£ be the distribution generated by the fields Di and 
d/dp^i, i = 1, . . . , n. The curves (|2.5p are now recognized as integral curves of a vector field 



X 



d ., ( d . r d 
dt +q 'd^ +Pa dp T a 



(2.10) 



which cancels the Pfaff system of differential forms ( "Cartan forms" ) 

^=dj£-K 4 d<f + i^di, (2.11) 

i.e., satisfies equations 

X_\u} r a = 0, for all r,a. (2.12) 

The system {w^} defines in J°° a distribution ("Cartan distribution" or CD), and the graphs T of solutions 
of (|2.1[) are the integral manifolds of this distribution, which is a necessary and sufficient condition for every 
surface in J°° to satisfy two requirements: First, to truly represent a graph of some function, i.e., to ensure 
agreement between coordinates p T a in J°° and the values of corresponding partial derivatives of a function 
p r (q_, t) represented by this graph. And second, to guarantee that this function satisfies Eq. (|2.ip [71IH1IS]- 

In the case of first order PDEs with m = 1, considered in a space of one-jets J\ : the trajectories, described 
by Hamilton equations, are the characteristic curves of a corresponding exterior differential system. These 
curves are uniquely defined in J\ , and so one and only one of them passes through every point on the graph T 
of the PDE solution. It is therefore natural to expect that the desired trajectories in J°° are the characteristic 
curves of the system This condition, however, happens to not be sufficient for the selection of q. Indeed, 

by direct calculation it is easy to check that the differential forms (|2.1ip satisfy the equation 

d^^dg i A^ J --f dtAut- (2-13) 
dpi 

From Eqs. (I2.13[) and (|2.12[) we immediately obtain that X_i &uj r a is a linear combination of forms (|2.1ip . 
Therefore, for every q the vector field X not only belongs to CD, but is also its characteristic field, and so in 
sharp contrast with the case of first order PDEs, there is a continuum of characteristic curves of CD passing 
through every point of T in J°° . 

It is instructive to consider the source of this difference. By analyticity, the graph of one and only one 
analytic function can pass in J°° through every point, which defines all spatial partial derivatives. Therefore, 
J°° is split into a foliation, each leaf of which is a graph T of the analytic solution of Eq. (|2.ip , so that every 
point of J°° belongs to one and only one leaf. These leaves are the integral manifolds of CD, so CD in the 
infinite jet space of analytic functions is completely integrableQ In contrast to J°°, in J{ CD is not completely 
integrable. Consequently, the dimension of the graphs T there is less than the dimension of CD, and so through 
every point of J\ different integral manifolds of CD (i.e. graphs of solutions T) pass, namely, with different 
values p\ k of the second derivatives on T. On the other hand, graphs T are formed by characteristic curves, 
so at every point these curves should belong to every graph that passes through this point. It turns out that 
this condition alone is sufficient to uniquely specify the characteristic field X at this point, including the value 
of q. Indeed, the corresponding PDE has the form 

%+H = Q. (2.14) 
at 



4 We note that this complete integrability cannot be considered a consequence of the Frobenius theorem and Eq. 112. 1311 . for the 
Frobenius theorem requires finite dimensionality of the space, and in infinite-dimensional space is no longer true [8] . 
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On a graph T of a given solution of (|2.14l) in J\ : the operator of the total derivative is 

D <=w +Pi h +P * k w k - (2 - 15) 

Similar to (|2.5p , the curves on T should satisfy a system of ODEs: 

p = Pl q l -H, (2.16) 
Pi = p\ k q k -D\H 

= &fa-^V(£+*£U (2.17) 



dpk ) V ®1 % dp 

Now, as characteristic curves should belong to every such graph, the dependence on p\ k for them must dis- 
appear. Recalling also that H does not depend on p, we then obtain from Eq. (|2.17j) the standard Hamilton 
equations for q % and pc 

■i dH . dH fn oN 

«=^' Pi = ~W (2 - 18) 

Thus, the fact that in J\ characteristics of CD have a Hamiltonian form is based on the specifics of J-J, 
namely, on the lack of complete integrability of CD there, which explains why this property cannot be gen- 
eralized to J°°. Fortunately, however, the requirement of being characteristics is not the only one which 
distinguishes Hamiltonian curves from all other curves that lie on the graphs of solutions of Eq. (12. 14)) in J\ . 
Let u> — dp — pidq l +Hdt be the (only) Cartan form in J\. The Hamiltonian curves in J\ may then be defined 
as integral curves of a vector field X, which cancels the 2-form duj, i.e., satisfies the condition X_i dco = [19] . 
This condition may of course be rewritten as £_l (X_i dui) = 0, for all £ G T^J-J). Now, in J°°, CD is defined 
by a sequence of forms {w£}, and a fruitful generalization of the above condition in J\ is to require 

£_l dul) = 0, for all r, a, £ g T£, (2.19) 

in J°°, which is a condition satisfied by the usual Hamiltonian curves of first-order PDEs with m = 1 after their 
prolongation from J\ to Jf°. In the following, we will call fields X that satisfy conditions (|2.19|) generalized 
Hamiltonian or, for short, simply Hamiltonian as they are called in J\. To find their form we observe that the 
internal product of X, given by the first line of Eq. (|2.10[) . with dui r a , is 

IJ ^ = E ~ E - ^ -Eif^ + EK-^ai*, (2.20) 

i v u rtri/ i StV u Pu SL , Pit 

where all summations are explicit, omits terms with (s, v) — (r, ai) for all i, and where we introduced 

functions V a = l r a (t : q, p, q) by the relations 

f^pUq'-K. (2.21) 

Now if the field £ belongs to T£, i.e. is a linear combination of Di and d/dp r ai , then only the first two terms 
in (|2.20[) contribute to £_l dcj£), and the condition (|2.19p gives 

? = JT^' ( 2 - 22 ) 
Pli = IU, (2-23) 
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where in Eq. (I2.22[) summation over r and a is not assumed. Equations (|2.23[) (with arbitrary a) just reproduce 
Eqs. (|2.5p (with a ^ 0), while Eq. (|2.22[) gives the desired expression for velocity^ For this expression to make 
sense, we must additionally require that the derivative there be independent of r and a. We will prove the 
following: dH^/dp r ai does not depend on a if H r satisfies the condition 

r) TJ r 

D k 1 — = 0, for all k,v ^ (2.24) 

(no summation over r and k here!), as, for example, if H r is linear in second and higher order derivatives of 
p r with constant coefficients. First we show that if any function F(t,q, {p}) depends on derivatives p r a only 
up to some finite order, and satisfies the condition D^F = 0, then it is a function of t and q J , j ^ k, only. 
Indeed, let r, v be such that dF/dp r a — for all a such that |<r| > \v\. Then applying the obvious identity 

■S- = Tr^-Dk - Dk-J^- , (2.25) 
dpi dp r wk dpl k 

where summation over k is not assumed, to the function F, we obtain immediately that dF jdp r v — also. 
Consequently, F cannot depend on any p r v at all, and is a function of t and q only, so that the equalities 
= DkF = dF/dq k prove our statement. Now for a function H r which satisfies condition (12.24[) . this means 
that for any nonempty multi-index v and any fc, dH r / dp r vk may only be a function of q J , j ^ k, and t. 
Consequently, H r has the form H r = H{ + where H{ = H{(t, q,p£, {p%, s ^ r}), and = ' s £j\ v \>\ a vV T y, 
where a v depends on t and q 3 , j £ v, only. It can then be easily seen that does not contain p r ai at all, while 
for any a ^ the only term with p r ai in H\ a is an additive term equal to p r ai dH[ j dp\ , and so all dFL r a j dp r ai 
are equal to dH\ldp\. 

We conclude that functions H r produce a Hamiltonian field if they satisfy the following "Hamiltonian 
conditions" (HC): First, for all r, the H r satisfy Eq. (|2.24|) (HC1). Second, given i, dH r jdp\ is independent 
of r (HC2). Equation (|2.22[) then gives for q % the values that are the same for all r and a, and, along with 
Eq. (|2.5p for p r al constitutes the desired system of equations of motion. We will now discuss the structure of 
the resulting theory in more detail. 

First note that the right hand sides of Eqs. (|2.22|) and (|2.5|) depend only on p r a with a ^ 0, and not on p r . 
Consequently, this system of equations splits into two subsystems. The first subsystem consists of Eq. (|2.22|1 
and Eq. (|2.5p with er ^ 0, which determine time histories (trajectories) q(t), p(t). The second subsystem, 
consisting of equations (|2.5j) with a = 0, shows how actions p r vary along these trajectories. The equations of 
the first subsystem, which is equivalent to the system (|2.22j) and (|2.23p , constitute the desired set of equations 
of motion. After the equations of motion are solved, equations of the second subsystem allow us to obtain 
expressions for the actions p r along corresponding trajectories by quadrature. As will be discussed later, if one 
makes equations of motion a starting point of a theory, then these expressions may be considered as defining 
actions for given trajectories. 

We will now show that Eq. (f2T23|) may be presented in a "Hamiltonian" form, similar to the second 



5 Note that for first order PDEs with one unknown function this expression does not depend on a and gives the same value of q' 
as in Eq. jSTDp . Consequently, one can bypass the geometric consideration above by simply postulating l|2.22ll in a general case. 
For all PDEs such that the right hand side of 112.2211 does not depend on r and cr, the resulting theory will be a generalization of 
the usual Hamiltonian formalism for first order PDEs with one unknown function. 
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Eq. (|2.18[) . Indeed, using Eq. (|2.22j) and making all summations explicit, we have 




(2.26) 



where v omits terms with (s, v) — (r, aj) for all j. Now consider H T a as a function of V r a ji 3 = • • • i n i 
and all its other arguments: H r a = H^(t, q, p',^^), where p' is the set of all pi such that v ^ and for any j, 
{s,v) ^ (r,aj). Let 9 be the graph of a ^-independent analytic function of q passing through a point with 
P-coordinates (t, q, p) in J°° . Let p' e (q) be the set of values of p' at a point of with a base coordinate q. 
Then using (|2.26p . Eq. (|2.23l) may be written as 

^ = -^^(t !q ,p' (q),py . (2.27) 

Although the right hand side of this equation is nothing but the last line of Eq. (|2.26p rewritten less explicitly, 
we prefer this form because of its obvious analogy with the standard Hamilton equation for momentum. 
Making the reference to graph implicit, the system of equations of motion (|2.22p and (|2.23|) can now be 
presented in the form 

g* = ^-fl;(t,q,p , (q)Xi), fi,-^:(i,q, P '(q),py. (2.28) 

We will informally refer to variables q l , p r a , and p r ai as forming an r-a sector of the theory. Evidently, in every 
r-a sector the theory looks like a standard theory of characteristics for the Hamilton- Jacobi equation with the 
Hamiltonian H^. in J\ — the same conclusion that may be drawn by comparison of Eqs. (|2.4p and (|2.14|) . or 
(|2.5p and (|2.16p . We will next explore other aspects of this similarity. 

2.4. Variational principles 

If one introduces "action forms" p — pi dq l — Hdtm J\ and p r a = p r ai dq l — dt in J°° , then the corresponding 
Cartan forms may be written as ui = dp — p and uf a = dp r a — p r a . For an arbitrary curve C in J-J, the difference 
of p at its ends is Ap — j c dp. Now if curve C lies on the graph F of the solution, then the tangent vector 
X at an arbitrary point of C cancels the 1-form to, and so X_i dp = X_i p. For such curves, therefore, the 
difference of p at their ends is given by the invariant Hilbert integral Ap — J c p [501 HI] ■ Similarly, in J°° , the 
difference of p r a at the ends of any curve C that lies on the graph of solution is 

Aj£ = I Pi- (2.29) 
Jc 

In J\ , the condition that the integral J c p remains stationary with respect to any variations of curve C that 
do not change the space and time coordinates of its ends may be used to select the curves (trajectories) which 
satisfy the Hamilton equations of motion [TO]. Similarly, in J°°, the trajectories that satisfy Eqs. (|2.28p may 
be selected by the condition that for every r and a, the integral J c p r a is stationary with respect to variations 
of C that do not change the base coordinates of its ends and are generated by any vector field which belongs 
to T£. Indeed, let C go from point A to point B in J°°, and let an infinitesimal e- variation transform C into 
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the curve C, which goes from A' to B' . By the condition stated above, the base coordinates of A' and B' are 
the same as of A and B, and so p r a — jS, p r a = 0, where the integrals are taken along straight segments 
connecting A with A' and B with B' . Consequently, the variation of the integral is equal to the integral over 
the closed loop AA 1 B' BA, which in turn, using Stokes' theorem, may be represented as an integral over an 
area D inside the loop: 







I Pi = 




lc 


J c JC 



Pi = f Pi = d P ;. (2.30) 

Jaa'B'ba Jd 

If C is an integral curve of a vector field X, parameterized by a base coordinate t, and the variation is generated 
by a vector field V, then the above integral is equal to e f t B V—l (A_l dp£)di. Since dp^ = — duj 7 a , the integrand 
here has the same form as the left hand side of Eq. (|2.19p . Consequently, repeating the derivation that follows 
(|2.19p we conclude that this integral vanishes for arbitrary field V 6 if and only if Eqs. (|2.22[) and (|2.23[) are 
satisfied (or, equivalently, Eq. (|2.28[) is satisfied). Clearly, this variational principle is completely equivalent to 
our original condition (|2.19l) . and may be considered to be its restatement. 

As in J-J, our trajectories may be also obtained from a different variational principle, a Lagrangian one. 
For that, we introduce the Legendre transformations of Hamiltonians H^, i.e. the functions 

dm 



dpi 



(2.31) 

*<" dp^ ° 

For Hamiltonians satisfying condition (|2.24l) . when a ^ 0, this transformation is of a trivial nature: As was 
discussed above, in these cases H r a contains p r ai only in an additive term which is linear in it, so in Eq. (|2.31l) 
this term gets canceled, and L r a is simply equal to the sum of the remaining terms with a minus sign. When 
cr = 0, however, H r does not have to be linear in p\. We assume that there are some values of r such that H r 
is a nonlinear function of pj, and will consider only those values of r below. For those values of r, velocities 
v l = dH r /dpi are nontrivial functions of momentums p[. These relations between velocities v % and momentums 
p\ are supposed to be resolved with respect to p-', expressing them through vasp[ = ip? (t, q, p', v) with some 
functions (p\, where, here and in what follows, p' is a set of all p s v such that v ^ and (s, v) ^ (r, i). The 
obtained expression for p\ should then be substituted into Eq. (|2.31[) with a — 0, resulting in the definition 

L r (t, q, p', v) = q, p', v) - H r (t, q, p', ^(t, q, p', v)) . (2.32) 

The above-defined Lagrangians L r can now be used to derive the principle of stationary action in its Lagrangian 
form from the invariant Hilbert integral. The derivation closely follows the one for the first order PDE [2"T] ; 
nevertheless, it is presented here for completeness and because of the complicating presence of (absent in the 
first order case) higher momentums. 

Let A and B be two points on the graph T of the solution that both belong to the same "true" (i.e., 
obtained by solution of Eqs. (|2.22[) and (|2.5[l ) trajectory Cab- Let q(t) be an arbitrary curve through the 
base, connecting the base projections of A and B. The difference between the values of p r in B and A may 
then be expressed, as in Eq. (|2.29l) . by the integral 

&Pab= I [^-ir(*,q,p',pD]di, (2.33) 

where C q is the image of the curve q(i) on the graph T. In this integral, p\ and H r may be expressed [55] 
through L r as 

# = ^ ? L'-(i ) q,p',v), H r (t,< l ,p',p r i )=p r i v i -L r (t,c l ,p',v), (2.34) 
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where the value of v % in these formulas should be set to a known value dH r (t, q, p', p\)/dp[, i.e., to the velocity 
of a true trajectory, passing through the corresponding point q. Substituting these expressions into Eq. (|2.33l) , 
obtain 

L r (t, q, p', v) + (<f - „') A L r {t> q , p ' ; v) 



dt. (2.35) 



Now the integrand here contains the first two terms of a Taylor series expansion of L r (t, q, p',q) in powers 
of q l — v i . When q % is close to v l , i.e., when the curve C q is close to the true trajectory Cab, we have 

f L r (t,q, J >',q)dt-Ap' AB = ± f (<f - v*)^ - V j )g^j L r (t, q, p', v)dt + O ((<f - v'f) , (2.36) 

which is of second order with respect to q — v and, therefore, with respect to the deviation of trajectory C q 
from the true one Cab- Consequently, the "action integral" J c L r (t, q, p', q)dt is stationary with respect 
to variations of C q around Cab, which is the way the principle of stationary action is formulated in J°°. 
The stationarity of the action integral is, therefore, an alternative condition, which may be used to select 
true trajectories, which satisfy Eqs. (|2.22[) and (|2.5|) . from arbitrary curves on T, which satisfy Eq. (|2.5j) only. 
Eq. (|2.27j) with a — then becomes the Euler-Lagrange equation in the usual way, and the difference between 
the values of the action function p r in two points may be expressed, as for first order PDE, as an integral from 
the Lagrangian along the true trajectory that connects these points: 



Ap r AB = / L r (i,q,p',v)dt. (2.37) 

JCab 

It is clear from Eqs. (|2.29[) and (|2.31j) that a similar expression may be written for /S.p r a with arbitrary r and c, 
the only difference being that when H r a is linear in p r ai , the corresponding Lagrangian L r a does not depend 
on v and so is a function of t, q, and {pi : v ^ 0, (s, v) ^ (r, ai)} only. Therefore, we always have 

&P r « AB = [ K^t. (2.38) 

J Cab 

2.5. Prom ODEs to PDE: infinite phase space formulation 

So far, we started with the PDE ([2~T]) and developed the system of ODEs (|2.28j) . Now we take system (|2.28|) 
as the starting point and will construct a corresponding PDE from it. The system lives in an infinite phase 
space V which has a geometry almost identical to that of J°°: The operator of total differentiation Di (|2.2[) . 
basis forms ui r a (|2.8j) . duality relations (|2.9j) . and Cartan forms uj r a (|2.11j) in V will be the same as in J°°, with 
the only difference that in V all multi-indices in these formulas and in all summations should be nonempty. 
Thus forms ui^ and ui r a exist only for a ^ 0, and integral manifolds of CD, defined by forms wj, are graphs 
of derivatives of analytic solutions of Eq. (|2.1j) . Also, vector field X, generated by Eq. (|2.28|) . is tangent to 
these graphs, has the same form (|2.10j) with q l = dH r jdp\ (no contribution with a = there!), and satisfies 
Eq. (|2.20|) . Equation (|2.38j) with a — then may be considered as defining the actions S A b, corresponding 
to arbitrary curve Cab, and for r such that H r is not linear in p\, Eq. (|2.27j) with a = will ensure that 
the Euler-Lagrange equations for these actions are satisfied, and so the principle of stationary action for them 
holds. The action functions S r (q, t) with given initial condition 5o(q) at t = to are defined as in classical 
mechanics [19] : Namely, let Sq = Oo-Sq, then for every space vector qo consider a trajectory in V that is a 
solution of Eq. (|2.28|) with initial conditions 

q(*o)=qo, P r M = S^(q ) , (2.39) 
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and define functions 5 r (q, t) by 

^(q,i) = SS(q )+ / L r dt, (2.40) 



Cab 

where Cab is the trajectory that ends at time t in a point B G V with a space coordinate q and qo is a space 
coordinate of a starting point A e V oi this trajectory. We will assume that trajectories don't intersect, and 
so this definition is unambiguous. While equations of motion (|2.28[) and their solutions describe individual 
trajectories, action functions S r (q, t) describe a family of trajectories selected by Eq. (|2.39[) . Consequently on 
trajectories that form the family, dependence of the initial momentums p^ = p^(to) of the trajectory on its 
initial space coordinate qo is given by 

p; o (qo) = S r CT (qo)- (2-41) 

For the just-defined functions S""(q, t), the following generalizations of classical results hold: First, for any 
cr^0 the functions S r a — d a S r satisfy 

S r M,t)=pl(q,t), (2.42) 

where p r a {<\, t) is a ^-coordinate of point B. Second, for any a, including a — 0, the satisfy the Hamilton- 
Jacobi-type equation (I2.4p 

J^(q,t) q,S) = 0, (2.43) 

where S is a set of all partial derivatives S£, a ^ 0, of functions S ,r . As in classical mechanics, Eq. f|2 .42[) means 
that the values p r a , which originally were independent variables evolving according to Eq. (|2.28[) . become also 
partial derivatives of action functions S r . 

As was the case with the principle of stationary action, the proof is similar to the standard one |19j , but 
we present it because there are some additional complications. Let Cab and Ca 1 b 1 be two close trajectories, 
with the base coordinates of A and A' being (qo,io) and (q ,to) and of B and B' being (q, t) and (q',t'). 
These trajectories are integral curves of the vector field X. Let action forms p r a for any a, including a — 0, be 
defined as in J°° . Now connect point A with A' and point B with B' by straight segments, and consider an 
integral of p r G along the closed loop AA'B'BA. By Stokes' theorem we have 



Pi = d^, (2.44) 

AA'B'BA JD 

where D is the region inside the loop. The difference V(t) between points of Cab and Ca'B' with the same t is 
given by the vector AA' , dragged (with parameter t— to) by a flow of the vector field X. Now, the vector AA' is 
the vector q — qo raised to the graph of the analytic function S'o(q) in V. Consequently, if q — qo = e l d/dq l , 
then A A' = e z Di, so A A' is a linear combination of It is easy to calculate that the Lie derivative of Di in 
the direction of X is 

[X,D i ] = -(D i q>)D j , (2.45) 

i.e., also a linear combination of Dj, and therefore, so is the difference V(t). Now, the integral on the right 
hand side of Eq. (pill) is equal to J t * B V(t)_l (A_l dp r a ) dt, but by Eq. ([2~2U)| X_i dp r a is a linear combination 
of Cartan forms that are canceled by any Di and, therefore, by V{t). We have, eventually, that for all a this 
integral vanishes, and with it the integral of along AA' B' BA, and so 

(2.46) 
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A B' 

But for a ^ 0, as in Eq. (|2.29l) . J B = p r a {A) - p r a {B), and similarly for J A , /£. On AA', t = t and so 
p r a = S r 0a , AA'_I eft = 0, and for all a 

A Pl= [ A plM = [ A SZ^dq* = S^(A') - SoM) , (2.47) 

I A J A J A 

which for o ^ is equal to p r a {A') — p r a {A). Consequently, for a ^ Eq. (|2.46p becomes 

/ S P^dg 4 - i^di = pl(B>) pl{B) , (2.48) 
which in the limit e % — s> 0, At — >• gives 

^|£(q,t)=J&(q,t), ^- t P r a (B) = -K(B). (2.49) 
Now if for all |<r| = AT, Eqs. (|2.42l) are true, then for such a and all i, 

SU(q,t) = -^S r M,t) = Apr (qt)=p r. (qt)) (2 50) 

and so f|2.42|) is also true for \a\ = N + l. Then for a = 0, since trajectories Cab and Ca'B' satisfy Eq. f|2.28j> . 
we have 

A pA 

p r = L r dt = S r (A) - S r (B) , (2.51) 

B JB 
B' a 1 

and similarly for J A , p r , while J A p r is given by Eq. (|2.47[) with a = 0. Eq. (|2.46p now gives 

/ p[dg l -i/ r di = S* r (S')-5 r (S), (2.52) 

JB 

or in the limit e 1 — > 0, A< — > 0, 

^5 r (q,t)=pJ(q»*), § i S r (B) = -H r (B), (2.53) 

which completes the proof of Eq. (|2.42l) . and then the second relations in Eqs. (|2.49l) and (12.53)) prove Eq. (|2.43p . 

We see again in Eq. (|2.40j) that solution of the PDE (I2.43[) may be obtained from solutions of the 
ODEs (|2.28[) . Conversely, any sequence of functions 5£(q, t) which satisfy 

-qT = SU, for all a, i, (2.54) 

and which also satisfy Eq. (|2.43p with initial conditions corresponding to a family of trajectories with given 
initial distribution of momentums (|2.4ip . may be used for integration of the equations of motion for trajectories 
of this family: Equation (|2.42p . read from right to left, gives for all t the distribution on the family's trajectories 
of momentums which satisfy Eq. (|2.28p . Indeed, we have for these momentums dp r ai /dt — didtS^. Using 
Eq. (|2.43l) . it is then easy to show that if a point q(f) moves with velocity q k given by the first equation 
in (I2.28p . then the time derivative 

#.(q(i),t) = ^ + ^rtg* (2.55) 



of a function p r ai at a point q(t) is given by the second equation in (|2.28[) . For Hamiltonians of first order, 
this is the basis of a Jacobi method of integration of equations of motion, and in the following we will call it 
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the "generalized Jacobi method" for arbitrary Hamiltonians. Now, when this method is used, it is obviously 
desirable to make it applicable to as large a class of trajectory families as possible. From this point of view, the 
formulation we used above is unnecessarily restrictive and may be generalized. Indeed, initial conditions that 
define the family's trajectories are given by Eq. (|2.41[) . In this equation, S 1 ^ are derivatives of the functions 
Sq. However, we saw that the sequence S G that is used in the generalized Jacobi method does not contain S r 
and includes only functions with a ^ 0. Therefore, all these functions are derivatives of «SJ, i — 1, . . . , n, 
while the functions S![ themselves and, consequently, their initial values Sq\, do not have to be derivatives of 
any other functions. On the other hand, we have from Eq. (|2.54[) that 



dqi 0y dq 



This means, that 1-forms S^dq 1 should be closed, d^S^dq 1 ) — 0, which will allow us to define the functions 
S r by 

S£(q)= r^d^ + ^Cqo), (2.57) 



and so S 7 Qi will be their derivatives. The integration in (12.57)) runs along arbitrary curves in configuration space 
Q that connect points qo and q, and the value So(qo), as well as the vector qo itself, are also arbitrary. Thus 
from the very beginning, the functions Sq are defined up to an arbitrary additive constant; moreover, they 
will be usual, single-valued functions only if configuration space Q is simply connected. If the fundamental 
group of Q is nontrivial, then in general Eq. (|2.57l) defines functions Sq as multi-valued, or single-valued on 
the universal covering space of Q. The branches of Sq may differ only by a constant, and so they all have 
the same derivatives Sq , a ^ 0. Therefore, these derivatives will be single- valued as they should be because 
the family has one, and only one, trajectory starting at every point of configuration space at t = to, and the 
functions Sq , a ^ 0, define initial momcntums of these trajectories. Consequently, in the currently considered 
statement of the problem, which starts with equations of motion in the infinite phase space V, the functions S r 
and their initial values Sq in the generalized Jacobi method are defined up to an additive constant, and in cases 
where the configuration space Q is not simply connected, may be multi-valued. Note that these conclusions 
are purely topological, not dynamical — they do not depend on the form of the Hamiltonians H T G or on their 
order. The simplest example is a family of trajectories on a circle that all have the same initial velocity v. The 
function Sq is then equal to vrip + const, where r is the radius of the circle, and ip is the angular coordinate 
on it. When d / 0, this function is multi-valued on the circle, but single-valued on the universal covering 
space R 1 . 

2.6. The case of complex-valued solutions 

We now allow complex- valued solutions of the PDE (|2.1|) . We only consider the case of one complex function 
p(q,t), the generalization to the situation when there are several of them being obvious. Let p 1 ' 2 be this 
function's real and imaginary parts, so that p = p 1 + ip 2 , and similarly H = H 1 + iH 2 . The conjugated values 
are p — p 1 — ip 2 and H = H 1 — iH 2 . By setting m = 2, the theory of the previous subsections may be applied 
directly to the functions p 1 and p 2 , treated as independent real functions with Hamiltonians H 1 and H 2 . 
However, it is often more convenient to express the same results via complex functions p and p, because in 
this representation they behave as if they were independent and also because the equations for p are simply 
the conjugated equations for p. 

As is usually done, introduce vector fields 

d I f d d \ d 1/9 d 



dp 2\dp l l dp 2 )' dp 2\dp 1+l dp 2 > 
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and 1-forms 



which satisfy duality relations 



dp = dp 1 + idp 2 , dp = dp 1 - idp 2 , (2.59) 



9 , d , 
— _l dp = 1, — — I dp = 0. 



9p dp 

— _ldp = 0, — _l dp = 1. 

op op 



(2.60) 



We have 



and so the operator of total differentiation may be written as 



a i 1 rot a i ro% a _ 

d<f op(j dp a 

Similarly, vector field X, Eq. (|2.10p . may be written as 



d d d 

°i = 1T7 + ^— + Poi ■ (2-62) 



x = l + « iDi - H °-k- S °lk- (2 ' 63) 



Since is an analytic function of pj , it satisfies the Cauchy-Riemann equations 

8H 1 _ dH 2 OH 1 _ dH 2 

9P] ' 'V; 9Pj 



(2.64) 



The first of these equations means that if p 1 ' 2 are considered as independent real functions with Hamiltonians 
H 1 and H 2 , then HC2 are automatically satisfied with the corresponding velocity 

1 = -TT = TT' 2 ' 65 

Expressing here d/dp 12 through d/dp and d/dp, and similarly H 12 through H and H, and taking into account 
that, being an analytic function, H depends on p 1,2 only through the combination p = p 1 +ip 2 , and H depends 
on p 1,2 only through p = p 1 — ip 2 , obtain 

dH 1 (d d\H + H lfdH dH\ .„ _ 

q M = W^WJ — = 2 [dp- + wj ■ (2 ' 66) 

It is easy to see that dH 2 /dp 2 gives the same expression for the velocity. 

For application to a theory of particles with spin, we also need to consider the case of a complex analytic 
function p(w, t) of complex coordinate w = w 1 + iw 2 with a Hamiltonian H{p w ), where p w — dp/dw. If w and 
p were real, the w's velocity would be equal to 

w=?—. (2.67) 

Op w 

It is easy to see that due to analyticity of all the functions involved and the corresponding Cauchy-Riemann 
equations, the same expression for w remains true in a complex case. Indeed, p w is given by the standard 
expressions 

Pw = pli + ip 2 w i = P 2 W 2 - ip X w 2 . (2.68) 
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From that, and using the Cauchy-Riemann equations, we have for the real and imaginary parts of (|2.67[) 

\d Pw J dp 1 ^ d P 2 wl ' \dp w ) d P 2 w2 dpl 2 

On the other hand, considering p 1,2 as independent real functions of real variables w 1,2 with Hamiltonians if 1 ' 2 , 
we have 

, ! _ dH 1 _ dff 2 . 2 _dH^ _dH^_ 
~ dpli ~ dp 2 wl ' W ~ d P l 2 ~ 8p 2 w2 ' 

where again we used the Cauchy-Riemann equations, and the first (resp. second) representation of p w in (|2.68|) 
for calculation of w 1 (resp. w 2 ). Thus, as it was for velocity <p , the HC2 for w 1 ' 2 are automatically satisfied 
due to the Cauchy-Riemann equations, and the lib 1 ' 2 are equal to 5R(9if / dp w ) and ^s(dH / dp w ) in Eq. (|2.69p . 
which proves (|2.67|) . 

2.7. Discussion 

The following general picture emerges from the above development. As for the standard case of first or- 
der equations, the PDE (12. 1|) with if r satisfying HC allows the introduction of a corresponding system of 
ODEs (|2.28|) . The values whose dynamical evolution is governed by this system are coordinates of a point, 
moving in a base, and partial derivatives at this point of unknown functions. On the other hand, as was just 
discussed, Eq. (|2.28|) may be considered on its own, as Hamilton equations are in classical mechanics, with 
the action functions being introduced later. In general, the system (|2.28[) is an infinite hierarchical system of 
coupled equations. We will not attempt its solution in this work; what will be important for us here is that 
this system has solutions whenever Eq. (|2.ip docs, and as for any system of first order ODEs, this solution is 
unique. The system then defines in J°° and V some trajectories, which lie in graphs T of solutions of (|2.1|) and 
are the characteristic curves of a corresponding exterior differential system. Like any characteristic curves [18] , 
these trajectories express solutions of Eq. (|2.1|) with given initial conditions and their derivatives as integrals 
of p r a or L r a dt along them. The direct proof of this statement, which doesn't use the theory of characteristics, 
is presented in section 2.5. The theory splits the whole jet space J°° into r-a sectors with coupled dynamics, 
described by Hamiltonians H^. The sectors share common coordinates q z , but there are no conflicts, because, 
thanks to HC, the dynamic they all define for these coordinates is the same. The structure of the theory in 
every r-a sector is similar to the one in J-J, but the value p r ai has a dual meaning: while on one hand, in an 
r-ai sector it plays the role of an action, obtainable from the above integrals, on the other hand, in sector 
r-a it is an i-th component of momentum, evolving according to the corresponding "Hamilton equation" with 
Hamiltonian H r a . In our equations, this duality may be seen especially clearly in the comparison of Eq. (|2.23j) . 
which describes the evolution of p r ai as an action in an r-ai sector, with the second equation of (12.28|) . where 
it evolves as the i-th component of momentum in sector r-a . As the second equation of (I2.28|) is just a differ- 
ent form of (|2.23p . sectors r-a and r-ai obviously agree on the dynamics of a variable p T ai which they share. 
A variable p r belongs to only one sector r-0, and so for its time derivative we have only the "action form" 
representation, given by Eq. ()2.5|) with a = 0. 

First order evolutionary PDEs with m = 1 satisfy HC automatically, and so the whole theory is completely 
applicable to them. They are different, however, from higher-order equations in the following important 
aspect: For any N > 1, the resulting equations of motion for q and p a , \a\ < N, form a closed subsystem, 
and the corresponding geometric theory may be formulated in a space of iV-jets Ji . Indeed, if the order of 
the PDE (|2.ip is equal to k, then H a contains the derivatives of orders up to \a\ + k. Since the Cartan forms 
w CT are expressed through p , which is of order \a\ + 1, and if CT , for first order equations, i.e., k = 1, and 



21 



any integer N > 1 the system of exterior equations {u a = , |er| < N} is closed with respect to the set of 
derivatives of p it includes. For this system, the space of iV-jets is sufficient, and the space Jf° is not 
necessary. On the contrary, for, say, second order equations, H a in the exterior equation w ff = contains 
the variables p v with \v\ = \a\ + 2. In order to ensure that these variables do indeed describe corresponding 
derivatives of the solution, which is represented by a graph in a jet space, we need to require that this graph 
also solves an exterior equation uj a > = with \<j'\ = \a\ + 1. But then H a i in w CT / will contain variables p V ' 
with \v'\ = \a'\ + 2 = \a\ + 3, and so the process will never stop, and the use of the infinite jet space Jf° 
becomes inevitable^ Similar considerations show that, while for a higher order PDE the expression for p a 
with (7/0 contains p v with \v\ > \a\, for a first order PDE it doesn't, and so for it a system of equations 
p a = l a with |it| < N is closed for any N > 1. The system corresponding to N = 1 is the simplest possible, 
but it still describes the evolution of the most important variables: q, p, and pi. A theory of this system in 
J{ is a usual Hamiltonian theory of first order evolutionary PDE, which is a part of their "full" , i.e. including 
all derivatives, theory, while the latter is a special case of our theory of satisfying HC evolutionary PDE of 
arbitrary order and with arbitrary m. 



3. Hamiltonian flow of quantum Hamilton-Jacobi equation 

Now we apply the technique developed above to non-relativistic quantum theory. We start with the one-particle 
case and consider the multi-particle situation later. By expressing the wave function as 

^(x,t) = ex P Qp(x,t)^ ) p(x,t)=£(x,t) + j#(x,t), (3.1) 

where p is complex and S and R are real functions of position and time, the one-particle Schrodinger equation 

in^=-^M + U{*,W (3.2) 
may be equivalently presented as an evolutionary PDE for p as 

| + tf = 0, (3.3) 
or as a system of evolutionary PDE for R and S as 

f + fl»-0, f + H"=0, (3.4) 
The Hamiltonian functions iJ, H s , and H R in the above equations are 

H - L$ +u+ ^L p ^ (3 - 5) 

H° = ^ + U-^(R] + R jj ), (3.6) 



H R = L( Sj R ] + l -S rj ), (8.7) 



where the indices denote corresponding partial derivatives and we extend the summation rule to expressions 
like p] =PjPj. 



6 Another reason to use an infinite jet space is Backlund's theorem l9l I10| , from which it follows that CD, defined by 1-forms 
1 12.1111 with Hamiltonians H r of higher-than-first order, cannot have characteristic fields in any finite jet space J k , k < oo. 



22 



Along with (ll.2p . Eq. p.3p and system p. 41) will be also called "quantum Hamilton- Jacobi equation(s)" 
(QHJE). Obviously, QHJE is equivalent to the Schrodinger equation (|3.2[) . and the action function p carries 
the same information as the wave function ip. In the following, for convenience, we will often discuss only 
one of these functions/equations, with the understanding that our conclusions may be applied, with proper 
modifications, to the other. Also, since there is only one wave function, we will always say "action function," 
even when there are several (two) of them. 

We can now see immediately that Eqs. p.3p and (|3.4p satisfy HC1. This is obvious in Cartesian coordinates 
used in Eqs. (|3.5I) - (|3.7[) . and is instructive to verify in the cylindrical and spherical coordinate systems. As 
they should (see section 2.6), Eqs. p.3[) and p.4[) also satisfy HC2 with corresponding velocity 



J — 



dH s dH R 1 



= -S jt (3.8) 



dSj dRj m 
which, in agreement with Eq. (I2.66p . may also be expressed as 

vi = ±-(p j +p j ). (3.9) 

Consequently, the theory of the previous section may be used. It means that in a space of analytic jets, 
corresponding to PDE p. 31) (or to a system of PDEs p. 41) ) . there exist trajectories, described by the system 
of ODEs (|2.28p . such that the solutions of the PDE and their derivatives may be obtained from the initial 
conditions by integrating the corresponding Lagrangians along these trajectories (see Eq. (|2.40p ) . The very 
existence of such ODEs and trajectories is just a mathematical fact, proven in the previous section. However, 
it raises an inevitable physical question: do the particles indeed move along these trajectories? Or, more 
practically: can peculiar features of quantum mechanics be understood, and its predictions reproduced, by 
assuming so? There are more questions. As we discussed, if the solution of the PDE is known, then the 
values of p r a and v may be obtained from it by the generalized Jacobi method. Therefore, the ODE and 
PDE formulations should be considered as two faces of the same theory, exactly like Hamilton equations and 
Hamilton- Jacobi equation in classical mechanics. But in classical mechanics, the roles of these equations are 
very different: while Hamilton (or Newton) equations provide the description of individual trajectories, the 
Hamilton- Jacobi equation describes an evolution of the (action) function which does not correspond to any 
particular trajectory, but is associated with a family of them. As was discussed at length in section 2.5, similar 
roles are played by the equations of motion (j2.28p and the PDE (|2.1I) in the mathematical theory of higher 
order equations. Now the other question is whether the situation in quantum mechanics is the same, so that 
the action or wave functions describe families or ensembles of trajectories, while the description of individual 
events/trajectories is provided by the system of ordinary differential equations of motion (I2.28p . In the rest 
of this work we defend a positive answer to these questions. As we already mentioned in the Introduction, 
we call this approach an analytical quantum dynamics in infinite phase space ("PAQD). We now start with its 
general description. 

As is clear from the previous section, the state of a particle at some moment t in "PAQD is defined by a 
triple (x, S,R) or, equivalently, (x, p,p). Here x is the position of a particle at time t, and S, R, p and p 
are sets of all the derivatives of the corresponding action functions in x at this time, so that S is a set of all 
S a (x,t) with ex 7^ and similarly for R, p and p. However, in the framework of PAQD, they are just a set 
of independent fundamental variables, identified by their multi-indices, which describe the state of a particle 
at time t exactly like components of momentum in classical mechanics. The evolution of a state is described 
by the equations of motion (|2.28p (the second of these equations is easier to use in the form (|2.23p ) , where 
the functions p r in that equation are now p 1 = S and p 2 = R or p 1 — p and p 2 — p. For future reference, 
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we present here expressions for Hamiltonians and some equations of motion. The first Eq. (|2.28p . i.e. the 
equation for x, takes the form of Eq. (|3.8[) for the (S, R) formulation and Eq. (13.91) for the (p,p) formulation. 
The operator of total differentiation for the (S, R) formulation is 

Dt = w +S «ik +R «ik- (3 - 10) 

Expressions for H s and H R are given in Eqs. f|3.6p and (|3.7p above. We also have 

1 Ti 1 ( 1 \ 

flf = — SjSji + Ui - — ( RjRji + - Rjji J , (3.11) 



H P - ~\ S i R J i + S J iR o + o $33* I ■ (3.12) 



— SjRji + SjiRj + — 
m\ 2 

The time derivatives of the actions S and R and their first momentums are 



S = ^-U + ^(Rl + Rjj ), (3.13) 
Si = -Ui + ^ + ±R jjl ), (3.14) 

* = -j^, (3-15) 



■R« — ( SjiRj + ^Sjji ) . (■'!..]()) 



For the (p,p) formulation, we only need equations for p a and H^, since the equations for p a and H a are 
obtained from them by conjugation in an obvious way. The Hamiltonian H is given by Eq. (|3.5p . and the 
operator of total differentiation by Eq. (|2.62|) . Then for Hi and Hik we have 

1 % . , 

* = m PjPji + i+ 2im Piji ' (3 ' 17) 

Sife = — (PjPjik + PjkPji) + Uik + -^7—Pjjik , (3.18) 

Tit Z/LTTXi 

while the time derivatives of the action and first momentums are 

p = i^—PjPj-U-^^—pjj, (3.19) 
2m 2vm 

1 h 
Pl = 2ro~ ^ Pj ~ Pi ' Pji ~ Ui ~ 2i^, Piji ' ^ 3 ' 20 ^ 

Pik = l^(pj-Pj)Pjik-^PjiPjk-Ui k -^pjjik. (3.21) 

It is not difficult to derive a general expression for H a . We say that a multi- index v is a subindex of the 
multi-index er, and write v C er, if there exists a multi-index /i such that er = z^/i. This multi-index [i will then 
be denoted as a \ v. Every multi-index is its own subindex, and the empty multi-index is a subindex of every 
multi-index. We also say that the multi-index v C a is chosen from the multi-index a if v is obtained from 
a in the following way: write a as a sequence of indices ij, ... then with this sequence fixed select \v\ 
members of the sequence to form v, and the others form a\v. Denote the summation over all such choices 
from a fixed sequence by X^xct" With this definition, it is easy to prove by induction that 

H ° = ^Y. V3»Pjo\v + u ° + ^ ' ( 3 - 22 ) 
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where the factor 1/2 accounts for the fact that in the sum over v every term appears twice. If not all indices in 
tr are different, then the same subindex v C a may be chosen from a in different ways. Consequently, there will 
be different choices that give the same (i.e., with the same v) contribution to the sum in (13.221) . For example, 
this will always happen when configuration space is one-dimensional, and the reader is encouraged to write 
formulas for H a with \a\ — 1,2,3,... in this case. It may be useful to present the summation in p. 221) in a 
form that contains only different contributions. Since the number of ways by which Vi indices i may be chosen 
from Oi of them in a multi-index a is equal to C^* = OiljvJ. (pi — Vi)\, the total number of ways by which a 
multi-index v C a may be chosen from a is C% = Yl7=i ^vj- Therefore, Eq. p.22j) may be rewritten as 

where the summation now is over all different subindices v of a. Correspondingly, the equations of motion for 
the p a , a 7^ 0, become 

Pa = i ^ ~ Vi)Pia ~ 2^ ^ C >JvP Jt y\v -U a - ^Pn* ■ (3-24) 

where the summation Ylvra excludes terms with v = and v = a . 

The emerging theory is in many respects similar to classical mechanics, but there are also important 
differences. As in classical mechanics, the particles in PAQD move along well-defined trajectories, with 
definite values of position, velocity, and all momentums at every moment of time. The states of the particle 
belong to an infinite phase space V, and the equations of motion (|2.28l) describe the evolution of these states in 
terms of Hamiltonian flow in P. For the (p,p) formulation, these equations take the form of (|3.9[) and (|3.24[) . 
and the corresponding Hamiltonian flow is generated by the vector field 



X 



d a d d j_ d 

t7 + v ~n + p° o — 1~ p& 

Cjt OX 1 OPa Op a 

d d d { ] 

= — + v i D l -H a - H a —, 

at dpc dp a 

where v\ p a , Di, and H a are given by Eqs. (|3.9[) . (|3. 241) . (|2.62p . and p.23p respectively, p a and H a are obtained 
by conjugation, and summation over a does not include a = 0. As in classical mechanics, initial value of the 
state uniquely determines its future evolution. The action function and QHJE are not needed for solution 
of the equations of motion. The action function is brought into use as an additional mathematical structure 
either by introducing a Taylor series (|2.6|) or via Eqs. (|2.39l) and (|2.40p of the previous section. For the (p,p) 
formulation, Eqs. (I2.40[) with initial condition p(xo,to) — Po( x o) take the form 

p(x,t) =p (x )+ f Ldt, (3.26) 

J Cab 

and conjugated equation for p(x, t), and similarly for the {S,R) formulation, they take the form 



S(x,t) - ^o(x )+ f L s dt, 

•'Cab 

R(x,t) = i? (x ) + / L R dt 

^ C A R 



(3.27) 



where L, L s , and L R are given by the right hand sides of Eqs. p,19[) . p,13p . and (|3.15p respectively, and Cab is 
the particle's trajectory, connecting points A = (xo,io) and B = (x, t). As was discussed in section 2.5, in the 
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spaces with a nontrivial fundamental group, the action function may be multi- valued, and it is always defined 
up to an additive constant. Consequently, only the derivatives of the action function are relevant, and so this 
function may be represented by a graph in P. It then describes a family of trajectories, determined by the 
given momentums at each position at some initial time. The same is true in classical mechanics; the important 
difference, however, is that while in PAQD Eqs. (12.39)) fix all momentums/derivatives, the corresponding 
classical equations [TH] fix only the first of them. As a result, in classical mechanics the action function cannot 
be considered as characterizing the individual state/trajectory of a particle: a given trajectory may belong to 
any of a continuum of different families of trajectories, with different action functions. Contrary to that, in 
PAQD, if the state of a particle belongs to some family, described by an action function, then by Eq. (|2.39p it 
determines all derivatives of this function at a point where the particle is. As an action function is analytic, 
it is equal, up to a constant, to the sum of a corresponding Taylor series. Consequently, in PAQD the state 
of a particle determines the action/wave function of a family, which includes it, and is, therefore, described 
or characterized by this function. This description, however, is not complete: since an analytic function can 
be expanded in a Taylor series at any point of space, there are different (i.e., with different q) members of a 
family that all have the same action/wave function. Thus a complete description of particle's state may be 
given either by a point in P or, equivalently, by a point in a base and an action function, defined up to an 
additive constant (or wave function, defined up to a constant factor). 

The action function introduced in this way satisfies QHJE, which expresses its time derivative through 
this function itself, regardless of which particular trajectory is responsible for its appearance, and so the 
action function obtains its own dynamics. Nevertheless, it is clear that in the framework of "PAQD, on the 
fundamental level of equations of motion, the action function is a useful, but purely mathematical entity: for 
determination of particle's trajectory, its use is neither necessary nor sufficient. However, the action function, 
or rather the wave function, gains physical significance when a family of trajectories described by it gains phy- 
sical significance. This will be the case when one considers the preparation of an experiment. Namely, as we 
will see later, using macroscopic control tools one can usually fix the wave function, but not the trajectory (i.e., 
not the specific "PAQD state) of a particle. This means, that with every macroscopically identical repetition 
of an experiment, the wave function of a prepared particle will be reproduced, but with a different specific 
trajectory. These trajectories belong to the just-described family, and make up an ensemble that the wave 
function is associated with. Thus the wave function reflects the preparation procedure and describes the 
properties of an emergent ensemble, but not individual events (trajectories) in it, in agreement with Einstein's 
views (see corresponding discussion in |23]). The statistical distribution of trajectories in this ensemble will 
be discussed later. 

When a particle moves in an infinite phase space, its position in configuration space moves with the velocity 
given by Eq. p. 81) . This is the same velocity that is attributed to the particle in the de Broglie - Bohm theory 
(DBBT) [TJjrrS], where the wave function and particle's position are considered as fundamental elements of 
physical reality. It is then postulated that the wave function evolves according to Schrodinger's equation (|3.2[) 
and guides the motion of a particle according to Eq. (|3.8p . Alternatively, it is assumed that relation (|3.8[) is 
satisfied at some initial moment of time, and then the particle moves according to Newton's law, but under 
the influence of an additional "quantum potential" , which is created by the wave function and is given by the 
part of H s , Eq. (13.61) . proportional to h 2 . With an additional assumption about initial statistical distribution 
of particles, DBBT is known to reproduce experimental predictions of QM. From the "PAQD point of view, the 
relation between PAQD, DBBT, and standard QM is as follows: While PAQD develops both the ODE part 
of the theory, describing the particle's motion, and the PDE part, which describes the evolution of the action 
function, the standard QM restricts itself to the PDE part, thus being an analog of the Hamilton- Jacobi part 
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of classical mechanics without its Newton/Hamilton ODE part. Consequently, to compensate for this missing 
part of the theory, QM employs the statistical interpretation, which postulates the missing part's results. The 
progress achieved by DBBT is based on the observation that the need for the interpretation disappears if one 
postulates just described dynamical law of particle's motion, for all experimental predictions of QM can be 
deduced from this law mathematically. However, in the absence of a full geometric picture and the theory of 
equations of motion, developed in section 2, this modification of the theory required a promotion of the wave 
function to the rank of a real physical field that guides the particle or acts on it (but is not acted upon) with 
a quantum potential. As a result, DBBT drew a picture of the world so alien to the generally accepted ideas 
about a possible structure of physical theory, that the majority of the physical community found it too hard 
to accept, in spite of the theory's success with some difficult issues of QM, such as the measurement problem. 
As was discussed above, far from declaring the wave function a real physical field, PAQD may deal without it 
at all. However, using the wave function may be convenient from the practical point of view. Thus for "PAQD, 
DBBT just implements the generalized Jacobi method: rather than solve the ordinary differential equations 
of motion, one can instead solve the Schrodinger equation or QHJE, and then get the particle's velocity from 
Eq. (|3.8[) , where momentum Sj is obtained from the real part of the action function by a simple differentiation. 
The same procedure works in classical mechanics, and so for "PAQD the DBBT program sounds exactly like 
a suggestion to consider classical mechanics as a theory of particles and real physical "action field" S that 
evolves according to the Hamilton- Jacobi equation and guides particles, forcing them to move with the velocity 
v = VS/m. Besides different physical picture, PAQD also differs from DBBT by an extra requirement of 
analyticity, which will become increasingly important in what follows. Nevertheless, the particles in PAQD 
move along the same "Bohmian trajectories" with velocity p.8[) as in DBBT, which will allow us to use, with 
proper modifications, some of its important results. 

The classical limit of PAQD is best seen in the (S, R) formulation. The Hamiltonian H s and the equations 
of motion for the action function S and its derivatives contain terms proportional to Ti 2 . When these terms 
are small compared to other, "classical" ones, they may be neglected. The equations for S and S a then 
decouple from the equations for R and R a , and directly turn into the system of equations for the theory 
with a Hamiltonian, given by the first two terms of H s , Eq. (|3.6p . This is a first-order Hamiltonian of 
classical mechanics, and the theory is classical mechanics, prolonged from the classical space of 1-jets J\ to the 
corresponding infinite jet space Jf° of the "full" theory, which describes all derivatives of the action function. 
As was discussed in section 2.7, in the infinite system of equations of this theory the standard equations of 
classical mechanics form a closed "classical" subsystem, which provides full information about the evolution 
of the action function S, its first derivatives, i.e., components of classical momentum, and, most importantly, 
the particle's position q. If the equations of the classical subsystem are solved, the higher derivatives S a can 
be obtained from the solution either by quadrature (|2.38[) or simply by direct differentiation of the action 
function S. Thus in a classical limit (or in a formal limit Ti — > 0) PAQD dramatically simplifies, both 
conceptually and in terms of its complexity, and reduces to this subsystem, i.e., to classical mechanics. 

4. The form of Hamiltonian, superposition principle, path 
integration, and wave-particle duality 

The evolution of the wave function ^(x, t) over an infinitesimal time interval e may be represented by a one-step 
Feynman integral as 

iPMH^ + 0(e 2 ), (4.1) 



exp 



1 t / x ~ y * 
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where L(v, y, i) = mv 2 /2 — U(y,t) is a classical Lagrangian, and A = y^2irihe/m is a normalization con- 
stant [53]. Let q(t) be some curve in the base space, and {p{t)} be the set of values of the action function p 
and its derivatives at time t at the point q(t): {p(i)} = {^(ft/i) ln^(q(t),tj }. At any time, we have then the 
wave function 



tp(x, t) = exp 



iE^(i)(x-q(t))' 



(4.2) 



At time f = 0, let the curve pass through a point q = with velocity v, so that q(0) = and q(e) = ve. 
Using Eq. (|4.2p and letting z = x ve, p a = p CT (0), and = p<r(e), we have from Eq. (|4.1[) : 



^^!^ zCT = 7 ln / exp 



i r / z — y \ j \ 1 



0(e 2 



(4.3) 



Schrodinger's equation is a consequence of Eq. (|4.1j) . therefore, Eq. (|3.3I) with the Hamiltonian func- 
tion (|3.5[) . and then Eq. (12.51) follow from it as well. It is, however, instructive to obtain that the evolution of 
momentums p a along the curve q(i) corresponds to Eq. (|2.5I) . i.e., that 

p' a =p a + e{p ai v l -H a ) + 0{e 2 ) (4.4) 



with H(j given by Eq. (|3.23l) . directly from Eq. (|4.3[) . For that, we need to find the coefficients of the expansion 
of the integral in (I4.3[) in powers of z 1 . It is convenient to introduce one more variable u = y — z and, using 
the expression for the Lagrangian, rewrite this integral as 



In 



exp 



im 
2he 



im 
2h 



V £ 



ie x - 1 



^ ^ cr! 



p CT (z + u) c 



n 



(4.5) 



The integral here is of the kind that may be evaluated using standard rules of the diagram technique |25j . The 
logarithm in front of the integral means that we should include only connected diagrams. The first term in 
the exponent defines a contraction (u 3 u k ) = —(eh/im) Sjk- Since we are only interested in the zero-order and 
first-order contributions of e, and the contraction is proportional to e, we have to consider only diagrams with 
one contraction or with no contractions at all. Then the expression in Eq. (I4.5[) will become the sum of the 
following contributions: The contraction of the term —{im/K)vcv in the exponent with itself gives —mv e/2, 
where 1/2 is a symmetry factor, and cancels the contribution of the third term in the exponent. As the 
potential term in the exponent already has a coefficient e in front of it, we can write there z CT instead of 
(z + u) CT , and then the contribution of this term to (I4.5[) will be equal to -eEatV 17 ') U a 7, a . In the last term 
in the exponent, it is sufficient to expand (z + u) CT up to the second power of u l . This term then becomes 
equal to (i/K) X) CT ( 2;Cr / '') (Pf +PaiU l + p G ijU l v? /2). Now the contribution to (|4.5[) of the term p a here is equal 
to J2a( z ' 7 / a *)Pcr: the contribution of the contraction of p a iU l with — (wn/7i)uv is equal to e ^^(z^/cr!)^^^ 1 , 
the contribution of the contraction of p„iU % with itself is equal to — (e/2m) J2vfj,( z " zfl I ^-l^VviPvAi where 1/2 
is a symmetry factor, and finally the contribution of the contraction of u 1 with u 3 in p a iju l u 3 is equal to 
— (efi/2im)^2 ty ('z' T /a\)pajj- Now, collecting all terms and comparing the coefficients for equal powers of z l in 
both parts of Eq. (|4~3]) . we obtain Eq. (|4~4"1) with H a given by Eq. (|3l23]> . 

We have thus demonstrated that the whole system of equations (|2.5[) with Hamiltonian (|3.5p may be 
compactly represented by one equation (|4.3[) . On the other hand, starting from the ODEs (|2.5[) with Hamil- 
tonian (|3.5[) , and reversing the above arguments, we can derive Eq. (|4.3[) in the framework of "PAQD. Equa- 
tion (|4.3[) . therefore, is equivalent to the system (|2.5p with Hamiltonian (|3.5[) . and, being augmented with 
the variational principles of section 2.4 for determination of the velocity v, may be taken as an alternative 



28 



starting point of the theory. Equation (|4.3[) then will fit in the general scheme of "PAQD as a separate 
postulate, restricting the possible form of the Hamiltonians in (|2.5p . For more general situations than the 
just-considered motion of a particle in a fiat space under the influence of a potential force, the form of the 
Lagrangian function in (|4.3[) will be different, for example, in magnetic field it will include a linear in velocity 
term (e/c)A- (x— y)/s, where A is a vector potential evaluated at the "midpoint" (x+y)/2 |26J. The Hamilton 
operator in Schrodinger's equation and the Hamiltonian function H in Eq. (|2.5[) will then be determined by 
this Lagrangian function in the same way as for the standard case above, and to be able to develop "PAQD 
we have to require that the function H satisfies HCf , Eq. (|2.24l) . In the development based on Eq. (I4.3p . this 
additional condition appears completely arbitrary and artificial. We will see in the next two sections, however, 
that it is also necessary for a derivation of the probabilistic interpretation of the wave function. 

As was just mentioned, the quantum Hamiltonian function in Eq. (|2.5[) . which describes the dynamics 
of a particle in PAQD, is determined by the Lagrangian function in (|4.3[) , or by the corresponding classical 
Hamiltonian, obtained from it in the limit e — > 0. This last Hamiltonian will necessarily be of first order, i.e., it 
will depend only on position and the usual momentums p.;, and not on any p a with \a\ > 1. Thus the approach 
that starts from Eq. (|4.3[) automatically reduces the variety of possible quantum Hamiltonians in (|2.5p , which 
were previously restricted by the Hamiltonian conditions only, to those which are obtainable in the described 
way from some classical Hamiltonian of the first order. As was discussed above, this classical Hamiltonian 
will then describe the classical limit of the corresponding quantum theory, which in turn will become its 
quantization. However, this quantization does not have to be unique. Indeed, while the quantum Hamiltonian 
in (|2.5p is determined by a Lagrangian function in (|4.3[) , which is written for finite e, the corresponding 
classical Hamiltonian is obtained from this function in the limit e — > 0. Consequently, there might be cases 
when different Lagrangian functions in (|4.3[) define different quantum Hamiltonians, but the same classical 
Hamiltonian in the e — > limit. For example, the quantum Hamiltonian in magnetic field would be different, if 
the vector potential in the term (e/c)A- (x — y)/e in Lagrangian was evaluated at other point than (x + y)/2; 
this ambiguity reflects operator ordering ambiguity in canonical quantization [26] . In such cases, we will 
regard these different Lagrangians as defining physically different quantum theories that nevertheless share 
a common classical limit. In other words, the more fundamental quantum theory must uniquely define its 
classical approximation, but not the other way around. In such situations, if competing theories are supposed 
to describe nature, then not more than one of them can do it right, and it should be chosen based on its 
phenomenological success. 

The next observation regarding the approach based on Eq. (|4.3[) as a foundation of the theory is that it 
automatically introduces the wave function, which in "PAQD is defined by Eq. (|4.2[) . as an object with a linear 
law of evolution (|4.f p . The superposition principle then follows immediately. As was discussed above, the 
particle's state in "PAQD may be described by its position and a wave function, and it is the wave function 
part of this description that is the subject of the superposition principle: if at some initial time t — to, the wave 
function ip, Eq. (|4.2p . is equal to a linear combination of other functions ip)-i k = 1 , . . . ,rik, of the form (|4.2p 
with the same q(to) as ip, then it continues to be that combination as time evolves. The position q(t) of 
the particle then evolves according to the equations of motion with the wave function ip. There is no such 
concept as superposition of a particle's position, and the time evolution of this position in a state with the 
wave function ip is not related in any simple way to evolutions in states with wave functions ipk- 

Equation (|4.3p is equivalent to Eq. (14.11) . from which Schrodinger's equation (|3.2p immediately follows. 
An even more important property of the approach that selects Hamiltonians in (|2.5I) using Eq. (|4.3p is that 
while the integral in f|4.3[) is mathematically well defined, and does not suffer from any difficulties that are 
usually associated with the path integration, the iteration of Eq. (|4.ip leads to a Feynman path-integral 
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representation of wave and action functions. Thus the value of the action function at some point of the base 
space may be obtained in two seemingly very different ways: either as action integrals p. 261) and p.27p along 
a particle's well-defined trajectory, coming to this point, or as a logarithm of a sum over paths. Obviously, 
this remarkable duality is a consequence of the fact that the quantum Hamiltonian (I3.5[) . which determines the 
particle's dynamics in VAQD, was obtained from the Lagrangian function of the path integral via Eq. (I4.3p . 
Still, it is not immediately clear how the action integrals p.26[) and Q3.27[) . which operate only with the values 
defined directly on a particle's trajectory, and not anywhere else, can reproduce the sum over paths, which is 
obviously affected by the whole neighborhood of the trajectory. The answer is that the trajectory itself, and 
therefore the action integrals, arc determined by an infinite system of (ordinary differential) equations, which 
depend on all derivatives p a of the action and U a of the potential. But the action function is analytic (this 
is one of the postulates of 'PAQD), and we also assume that the potential function is analytic (and believe it 
always is analytic in nature). Consequently the trajectory, using these derivatives, obtains the full knowledge 
of the action and potential functions everywhere, and with it the ability to reproduce results obtained by path 
integration. In other words, the action integrals p.26[) and (|3.27[) utilize the information about analytic action 
and potential functions which is contained in their derivatives at the points of the actual particle's trajectory, 
while summation over paths uses the values of these functions on the whole space directly. But path integration 
provides a purely wave description of a particle's behavior, which naturally explains such characteristically 
wave phenomena as interference and diffraction. Therefore, it is because of a special form of the Hamiltonian, 
obtained from Eq. (|4.3p . and analyticity of the action/wave function, that the particle, which moves along 
a single trajectory, exhibits at the same time the characteristics of a wave, thus possessing the property of 
wave-particle duality. For example, in agreement with conclusions of [2 7) , in a two-slit experiment the motion 
of a particle, passing through one slit, may depend crucially on whether the other slit is open or closed, even 
when the difference between the classical forces acting on the particle in these two cases is negligible. This 
is of course a purely quantum effect, completely impossible in classical theory, where a particle's trajectory is 
determined by a finite system of equations that depend only on the first derivative of the potential. Therefore, 
the wave-particle duality in the quantum domain receives a simple and natural mathematical explanation 
in PAQD. 



5. Variational principle, continuity equation, and invariant 
measure 

Equations (|3. 41) . p.7p for R(x,t) may be rewritten in the form of a continuity equation 

-|- + divj = 0, (5.1) 

where j° — \ip\ 2 — e 2R , j = j°v, and v = VS/m is the particle velocity. The invariance of the measure, 
associated with conserved current (j°,j) = |V>| 2 (l,v), is used in the next section to demonstrate that the 
particle's probability density is equal to \ip\ . Therefore, the conservation of this current is a very important 
element of the theory, and in this section we present several different proofs of it, which will allow us to better 
elucidate its origin. We will also derive an important expression for a corresponding invariant measure. 
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5.1. Variational principle and current conservation 



Equation (|4.ip implies that Schrodinger's equation may be obtained from a stationary action principle. Indeed, 
consider the value 



H(^*) 



de 



■0*(x) exp 



■eL 



x-y 



y,t 



dx % dy % 



(5.2) 



Since the integration with a factor exp(ieL/fi) propagates ip(y) from time t to time t + e, and ifi*(x) from t to 
t — e, we obtain, obviously, 



V(x) 



■H(^,V*), 



(5.3) 



<5^*(x) v " w 5V(x) 
so that H(?/>,V'*) is a Hamiltonian function with respect to the canonical field coordinates ?A( X ) an d their 
conjugate momentums ^>*(x) ( or coordinates V'*( x ) an( i momentums — V>( x ), which differs just by a canonical 
transformation). Being Hamilton equations, Eqs. (j5.3[) , which are equivalent to Schrodinger's equation and 
its conjugate, follow, after standard discretization, from a stationary action principle in a Hamiltonian form 
6jL(if},ip*)dt = 0, where 

h 



L(V,V* 



v>*(x)^( X );q^-h(^*) 



(5.4) 



(see [19] and section 2.4) and the factor Uji is introduced for convenience. To calculate H(ip,ip*), integrate 
over JX m Eq. (|5.2|) to get 



where 7J is a Hamilton operator, and so L(^, V'*) = J^-Ili where the Lagrangian density £ is 



£ = ^*(x) ^- d t + H)^(x), 
For the standard Hamiltonian of Schrodinger's equation (|3.2j) we have then 



£ = V*( X ) - — dj + U U(x) 



(5.5) 



(5.6) 



(5.7) 



i 2m 3 

As is well known [28], both equations of motion and conservation laws, considered below, remain invariant 
with respect to adding a total divergence to the Lagrangian density. By adding appropriate terms to the above 
expression, we then obtain a familiar variational principle 6 J L J\ i dx l dt = with symmetric Lagrangian 
density 

c = - - ^d t r) + 7T- 9ii>*d^ + un> , (5.8) 

2% 2m 
and it is indeed easy to verify directly that the corresponding Euler-Lagrange equation 

dC 



dC 



i=0 



0. 



(5.9) 



is equivalent to the Schrodinger equation 

By Noether's theorem [7J[TUJ[3H], if a transformation ip —> ip + aA, ip* ip* + aA* with infinitesimal real 
parameter a changes the Lagrangian density C just by adding a total divergence to it, 



(5.10) 
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then the solutions of the equations of motion (i.e., in our case Schrodinger's equation Q3.2p ) satisfy a local 
conservation law 

d o j° + d kJ k =0 (5.11) 
with a current j l which is, for a first order Lagrangian C, equal to 

As was discussed at the end of section 2.5 and in section 3, by their very construction the action function 
is defined up to an additive constant and the wave function up to a constant factor. Therefore, we should 
expect that the corresponding transformation does not change the equations of motion, and so the original 
and transformed Lagrangian densities differ by a total divergence only. Since the wave function is complex, we 
should consider two different transformations. Under the scale transformation if> — > e a tp, tp* — > e a ip*, when 
the wave function satisfies the Schrodinger equation, the Lagrangian density C indeed changes as in (j5,10[) 
with 

(A°,A)= (o, Aty* V ^ + W)Y (5.13) 

but the sum of the first two terms in Eq. (|5.12l) in this case is equal to A 1 , and so the total current f vanishes 
and the scale invariance does not lead to any conservation law. The phase transformation ip — > e la / h ijj, 
is more useful. Lagrangian density C is invariant with respect to it, i.e., satisfies Eq. (I5.10P 
with A* = 0. Consequently, the corresponding current 

(i°,j) = U*1>, ^ {rvi) ~ (5 - 14) 

is conserved, i.e., satisfies Eq. (|5.11[) . which in this case coincides with a continuity equation (|5.1|) . Therefore, 
the current (j°,j) is equal to |^| 2 (l,v), where the velocity v is given by Eqs. p.8[) or (|3.9p . To obtain this 
form of the conserved current, we used an explicit form of a Hamiltonian here. However, in section 5.2 we 
will show that this result has a much more general character, namely, for a wide range of possible Lagrangian 
functions L in Eq. (|5.2p . the current, which is conserved due to the phase invariance, is equal to |^| (l,v) 
with v given by Eq. (|2.66p (which for the standard Hamiltonian coincides with (|3.9p ) . 

5.2. Current conservation for Hamiltonian operators of general form 

It is desirable to derive the conservation of the current (j , j) = |"0| 2 (1, v) under more general assumptions than 
above where we used an explicit form of a standard Hamiltonian. Here we will show that this conservation 
follows from the phase invariance of the Lagrangian density C, Eq. (|5.6p . for an arbitrary quadratic in velocity 
Lagrangian function L in Eq. (I5.2p . provided it satisfies some simple conditions. 
First substitute into (|5 . 6[) the representation (|3.ip to get 

£= \iJj(x)\ 2 ( Pt + H) , (5.15) 

where the Hamiltonian H, which corresponds to the Lagrangian function L in (|5.2p . is a function of the space 
derivatives of the action function p. We have from (|3.ip and (|5.2I) 



as 



6=0 



W — — In / cxj) 



H £L+ H P{X 



ndu 1 
(5-16) 



where we introduced u = y — x. As in section 4, expand the exponent in (|5.16p in powers of u l and consider 
W as a generating function for connected diagrams. We will assume that similar to the case of a standard 
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Hamiltonian, the corresponding contraction (v?u ) is purely imaginary and proportional to e. To account for 
a possible presence of a magnetic field, we allow the product eh to have a vector potential term — (e/c)AkU k 
inside it [55], but assume that there are no other e- independent and linear in u k terms there. The Lagrangian 
density (|5.6p will then be a function of the wave functions ip and ip* and derivatives of ip up to a second 
order. Using the corresponding formulas for second-order Lagrangian functions [71 1101 [55] . the current, which 
conserves due to invariance of the Lagrangian density C, Eq. (|5.6[) , with respect to the phase transformation 



ip -)• e la/h ip, ip* -)■ e" 



v//i 



will then be equal to (j°,j), where 
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(5.17) 



h \dip k dip k i, 

We now want to rewrite these expressions in terms of derivatives of H over pk and pki ■ 
We have from Eq. (|3.ip the following relations between partial derivatives of p and ip: 

h ip k 
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(5.18) 



Let pki and pik enter the expression for C symmetrically. We have then from (|5 . 1 8|) 

dC hi ( dC 
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(5.19) 



dip, 



where the factor of 2 in the first equation compensates for the dropped contribution of dC/dpik- Substituting 
these expressions into Eq (|5.17l) and using Eq. (|5.15[) and condition (|2.24[) . we obtain for the current 
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(5.20) 



h dpki . 

so that j° has the right form, and we need to evaluate the derivatives of H in an expression for j k . For every 
function B of coordinates u l , we denote by (B) the corresponding sum of connected diagrams produced by 
the generating function W: 



(B) = 



exp 
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p + p k u k + ^Pkiu k u l 



du l 



exp 
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- eh + - [ p + p k u K + - p k iu k u l + 



ndu 1 



(5.21) 



where the derivatives p a are taken at the point x where the Lagrangian density (|5.15|) is evaluated. We have 
then from Eq. (|5.16j) 
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e=0 



where we used the fact that the only nonzero contribution to d(u k )/de for e = comes from contraction of u k 
with e-independent terms in the exponent. On the other hand, we have 
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and so 

dH „ i / e . \ 

jr- = 2- W - -At — , 

(5.24) 

— = 2- [pi--Ai) — , 

where to obtain the second equation we conjugated the first one and used that dH/dpki is purely imaginary. 
Now using these equalities in Eq. (|5.20[) for j h , we obtain 

where v is the particle's velocity (I2.66|) . as was required. 
5.3. Invariant measure 

For an arbitrary current (i° , i) in the base space M that satisfies a continuity equation dti° + dki k = 0, a form 
v = i°ft, SI = dx 1 A ■ ■ • A dx™, integrated over any subspace of M corresponding to a fixed time t, is invariant 
with respect to a vector field Y — dt + u k dk with u = so that this form defines on M a measure that is 
invariant with respect to the flow of Y. Indeed, let D be an arbitrary cell in configuration space, every point of 
which moves with velocity u. If ds is an element of the boundary of D, orthogonal to a unit vector n pointing 
outside, then by the continuity equation over a time interval dt the integral J D v will reduce by indsdt due 
to the current i through ds. On the other hand, since the element ds moves with velocity u, during time dt a 
volume undsdt will be added to D, and with it a value i a undsdt added to J D v. Therefore, this integral will 
not change, and so v is invariant with respect to Y. 

The same result may also be obtained by direct calculation. Indeed, using equalities Y(y) = d(F_l v) + 
Y_i dv and dil = 0, we obtain 



(5.26) 



Y(y) = d(i°r_i ci) + y_i d(i°n) 

= di° A (Y_i CT) + i°d(Y—l Q) + Y—l di° A fl . 
In the last term, substitute F_l di° A fi = Y(i°) fl - di° A (Y_\ ft) to get 

Y{v) =i°d(Y_ifl) + Y{i°)n. (5.27) 
In the second term of this equation, we have Y(i°) = dfi° + u k dki°, and in the first term 

Y _j £1 = ^(-l) fe - Vtda; 1 A ■ • • A dx™)', (5.28) 

where ( )' means that the factor dx k in the product is dropped. From this, we have 

d(Y_l Q) = n d k u k + dtJ2(-l) k ~ 1 d t u k (dx 1 A • • • A dx")'. (5.29) 

k 

We will integrate Y{v) over surfaces with fixed t in the base space M, so we are interested in a pullback 
7T( Y(f), where n t : Q — > M is the natural embedding that maps configuration space Q into such surfaces. 
Then obviously the term with dt does not contribute to such integrals, 7r t *di = 0, and collecting remaining 
terms we have n% Y{v) = il (d t i° + u k dki° + i°dkU k ). But the last two terms sum to dk{u k i°) = dki k , and so 
from the continuity equation we obtain tt^Y(i/) = 0. 

The vector flow X in V AQD is defined by the first line of Eq. (|3.25p . Its first two terms have the form of 
the just-considered vector field Y with respect to a current (|5.14[) , and so conserve the form u — j°Cl, while 
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the last two terms, where a ^ 0, when acting on this form give zero. However, the flow X is defined in the 
infinite phase space P, rather than in the base space M . Consequently, to formulate the invariance condition, 
we use the map ir t : Q — > P, which projects configuration space Q into a part of the graph of a solution of 
Schrodinger's equation with given t. With so-defined 7r t , we have then the desired identity tt^X(oj) = 0, which 
expresses the invariance of the form uj and corresponding measure with respect to the Hamiltonian flow (|3.25p . 

6. Probability density 

In the previous sections we studied the mathematical structure of PAQD. Here we start considering its physical 
implications, i.e., experimental consequences of the assumption that particles move along the trajectories that 
we discussed. It is then natural to think, and we will confirm it later, that in PAQD the experimentally 
measured particle's position should be equal to its position in P. In this case, what can PAQD say about the 
distribution of this measurement's results? 

We believe that in every repetition of an experiment, in which the particle is described by a wave function "0, 
its position coordinate in P assumes a random value, determined by a specific history of this particular 
repetition. In the mathematical limit of an infinite number of such repetitions, the results form an ensemble 
that determines the particle's probability density p: the probability of finding the particle in any volume 
element is equal to the relative number of ensemble members with the particle inside that element. It seems 
natural to assume, in agreement with experiment, that this probability density is determined by the wave 
function only, and not by the way in which the ensemble with this wave function was created. Then once 
created at some time t, the ensemble remains representative for all future time, for one of the ways to create 
an ensemble at any time t' > t is to create it at time t and let it evolve till time t' . But during such 
evolution, every volume element dV, propagating with 'PAQD Hamiltonian flow, continues to contain the 
same ensemble members, and so the probability pdV to find the particle in this element remains constant, 
which means that the change in the probability density p in the element is inversely proportional to dV. On 
the other hand, as was demonstrated in the previous section, the product jf/') 2 ^ in this element also remains 
constant. Consequently, along any given trajectory, the probability density p should be proportional to \ip\ 2 - 
The coefficient of proportionality can, by this reasoning, depend on trajectory. However, we note that once 
selected, this coefficient should remain fixed in the presence of any external fields that may be applied to the 
particle in the future. Since such fields can shuffle trajectories in an arbitrary way, but the coefficient should 
remain a continuous function of trajectory, it is clear that for all trajectories it must be the same. Moreover, 
even if for some reason configuration space Q splits into two subspaces Q\ and Q2 such that trajectories never 
cross from one of them to the other, according to the way the wave function is brought into 'PAQD, that will 
only mean that rather than being defined up to one constant factor in Q, the wave function is now defined 
up to two independent constant factors in Qi and Qi. Obviously, these factors can be chosen in such a way 
as to make the coefficients of proportionality between p and |^| 2 in Q\ and Q2 equal. We conclude that the 
probability density p should be proportional to \ip\ 2 with a constant coefficient, or equal to it if J \ip\ 2 dV = 1, 
i.e., if the wave function is normalized. Thus in PAQD this classic relation between the wave function and 
probability density becomes just a property of the wave function, rather than its main physical meaning. 
It is worth recalling that this property is a consequence of such fundamental elements of the theory as the 
possibility of obtaining the Schrodinger equation from a variational principle (which in turn follows from the 
quantization via one-step Feynman integral — see the corresponding discussion at the beginning of section 5) 
and the definition of an action function, which leaves the freedom to add an arbitrary constant to it and which, 
therefore, for the theory following from a variational principle, results in a corresponding current conservation. 
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Also, in the process of proving that the space part of the current has the desired form (|5.25[) for a general 
Hamiltonian function H obtained from the one-step Feynman integral, we had to require that this function 
H satisfies HC1 (I2.24[) . Thus for any theory based on the one-step Feynman integral, the condition (|2.24[) is 
needed for the probabilistic interpretation of the wave function as well as for the very possibility to develop 
■PAQD in the first place. 

The same expression for the probability density may be also obtained in a different manner, if we count 
the number of possible ways by which a given ensemble can be created. This means the following: We break 
configuration space Q into small cells, so that in every cell the probability density p can be considered constant. 
The ensemble of N points, representing the particle in a state with a wave function tp, is described by the 
numbers Ni of points in every cell, Ni — N. The series of N experiments, which form an ensemble, is then 
characterized by the sequence i\,...,ijsi of cells the particle was found in in each experiment, so that such 
ensemble can be created in Af = N\/Y[ i Ni\ different ways. These numbers Ni and consequently A/", depend 
on the specific way the space Q is split into cells, and so some reasonable prescription for a way the splitting- 
is done should be made. We will demand, as we did above, that once created, the ensemble should remain 
representative for all future time. The splitting of configuration space Q into cells, therefore, should be such 
that once it is done and fixed, the number J\f for every ensemble, i.e., the number of ways this ensemble can 
be created, remains constant with time. But the image of every cell i, corresponding to the particle's flow over 
arbitrary time t, has the same invariant measure |^>| 2 dV and contains the same Ni ensemble members as the 
cell itself. It is then clear that the demand will be satisfied, if in the limit of vanishing cell volumes, they will 
all have the same invariant measure, which we denote as Ar. For our ensemble, described by the distribution 
of points in configuration space, we need now a characteristic of the number of its realizations that remains 
finite as N — > oo and Ar — > (in this order). The number of cells in a small area of configuration space Q of a 
particle is proportional to 1/ AT, and in a small area of configuration space Q N of the ensemble — to l/(AT) N . 
As Af is proportional to this number we need to factor it out, i.e., to consider (AT) N Af. The function of this 
last number that stays finite as N — > oo is the Nth root of it, and so we come to considering ATAf 1/N . The 
value 

S G (p) = lim lim InfArM 1 ^) 

/ ! X (6.1) 

= lim lnAr+ lim — In AT 



Ar^O \ Af^oo N 

will be called the Gibbs entropy of the ensemble, representing probability density p in configuration space Q. 
We have then, using Stirling's formula, 

lim -^\nAf = lim j- In , = - lim V ^ In ^ . (6.2) 

N^oo N N-yoo N U NA JV-s-oo-^ N N 

i 

Further, limjv-^oo N/N = p i AV l = {p l /\%p i \ 2 )AY and Pi^Vi = lim^oo E l N i/ N = h therefore 



S G (p) = A lim o 



InAr-^A^lnj^Ar 



AVi 



lim YAViPihx^- (6.3) 
v.^o ^ pi 



p 

In the absence of circumstances that make some cells preferable compared to others, every sequence i±, . . . , 
should be assigned equal probability. In the limit N — ► oo, the emerging ensemble ("Gibbs ensemble") should 
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with < is 



M 2 

lim —r = lim exp 



then maximize TV or Sg(p)- Indeed, in the limit N — > oo, the relative frequency of emergence of two ensembles 

^(4 2) "4 1) )] =0- (6.4) 
Now, since Inx < x — 1, we have from Eq. (|6.3[) , when the wave function is normalized, 

Sa(p) < J p(^--l\dV = 0, (6.5) 

Sc(p) achieving its maximum possible value of for p = |i/>| 2 , which will, therefore, be the observed probability 
density. 

As was discussed in section 3, the particles in PAQD move with the same velocity (|3.8|) as in the theory 
of de Broglie - Bohm. It was shown by Bohm [12] that this law of motion preserves the standard quantum 
form of the probability density: if p is equal to |^| 2 at some initial time to, then it will stay equal to it for 
all t > to. It was hypothesized p~21 [14] that an arbitrary initial distribution would converge to the stable 
density \ip\ 2 for t of the order of some "relaxation time," in the same way as macroscopic systems converge to 
thermal equilibrium. The derivation, presented above, shows that this hypothesis is unnecessary. Nevertheless, 
especially because we are using the concept of entropy, and looking for a distribution which maximizes it, it 
is instructive to compare the situation in PAQD with that in classical statistics. We present a brief sketch of 
statistical distribution and entropy growth in classical statistics, based mostly on the works 29, 30 , in a form 
convenient for such comparison in the Appendix. From the discussion there, the following conclusions may be 
drawn: 

- The convergence to thermal equilibrium in classical statistics is related to the growth of Boltzmann en- 
tropy Sb, rather than Gibbs entropy Sq, which is maximized in 'PAQD. 

- The growth of Boltzmann entropy is related to such properties of macroscopic systems as possibility of their 
crude, but adequate, description; as typicality (i.e., practical equality of observable magnitudes of additive 
physical values to their averages over microcanonical ensembles); and as possibility of replacement of one 
ensemble by the other in the process of these systems' time evolution (see details in the Appendix). These 
properties exist only in macroscopic systems that consist of enormous number of particles, and don't have any 
analogs in one-particle dynamics, classical or quantum. 

- The nature of the quantum distribution \ip\ 2 is identical to that of the microcanonical distribution in classical 
statistics. Both distributions maximize the corresponding Gibbs entropies, and emerge not because of the large 
number of particles in a system, but because of the infinite number of systems, be they one- or multi-particle, 
in the Gibbs ensemble. According to Eq. (|6.4|) . the ensemble with less than maximum Gibbs entropy has 
zero probability to arise. Consequently, all observed distributions automatically have the maximum possible 
values of their Gibbs entropy, in contrast to macroscopic systems' Boltzmann entropy, which grows due to the 
physical process of thermalization. 



7. Multiparticle systems and quantum particles in a macroscopic 
classical environment 

We now extend our approach to multiparticle systems. For a system of n p particles in n-dimensional space, 
we do it by directly combining n p one-particle n-dimensional configuration spaces into n p n-dimensional con- 
figuration space Q of a system. The theory of previous sections will be generalized in a straightforward way to 
look like a one-particle theory with corresponding Hamiltonian in an n p n-dimensional space Q. In particular, 
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the wave function of a system is related to the corresponding action function as in Eq. (|3.ip : 



i 



V>(xi,...,x np ,i) = exp I- p(xi,...,x„ p ,i)J , 
p(xi, . . . ,x„ p ,t) = 5(xi, . . . ,x„ p ,i) + j i?(xi, . . . ,x„ p ,i) , 



(7.1) 



and momentums in the infinite phase space are partial derivatives of the action function with respect to the 
components of Xi , . . . , x„ p . The Hamiltonian function (the multiparticle analog of the one-particle Hamilto- 
nian ([3~5]) ) 

n p 2 * n p 

H = yi^ + u+-y^ : (7.2) 

^ 2m k i ^ 2m k v ' 

fe=l K fc=i K 

where summation over repeating indices jk is from 1 to n, is obtained from the multiparticle Schrodinger 
equation in the same way as in section 3, and defines the evolution of momentums and particle velocities by the 
equations of motion (|2.28j) . so that in particular the velocity of the fc-th particle is = Vfc5(xi, . . . , x„ p , £)jm k . 

For macroscopic systems, the part of the action function related to their directly observable macroscopic 
degrees of freedom is much larger than Planck's constant Ti. As was discussed at the end of section 3, in the 
corresponding equations of motion the terms with h may be dropped, and then these equations reduce to those 
of classical mechanics, so that these degrees of freedom will exhibit a classical behavior. The wave function 
describing these classical degrees of freedom is given by Eq. (|7.1[) with the action function S that solves the 
multiparticle analog of the classical part of Eqs. (|3.4[) . (|3.6I) . i.e., the classical Hamilton- Jacobi equation 



where the macroscopic degrees of freedom are combined into the vector q, H c (q, p) is a corresponding classical 
Hamiltonian, and where by derivative with respect to a vector we understand a vector made from derivatives 
over the corresponding components. As was discussed in section 5, the wave function amplitude A = e R always 
satisfies a continuity equation, which in this case has the form 



OA 2 A 



A 



2 9H c (q, p) 



dpi 



p—dS/dq 



(7.4) 



The standard quantum- mechanical derivation of equations (|7.3I) . (|7.4p for the action function and amplitude 
in the quasiclassical case may be found, for example, in |31j . 

We now want to consider a combined system, consisting of macroscopic objects interacting with quantum 
particles. The same consideration applies to the interaction of macroscopic objects with their own internal 
(like electrons' or phonons') microscopic degrees of freedom, which should be described quantum-mechanically. 
In fact, it will be sufficient for our analysis to consider an extremely simplified situation where a macroscopic 
object is represented by one particle with a large (macroscopic) mass M in the limit M — > oo interacting with 
a quantum particle with a fixed (microscopic) mass m. Let the Hamiltonian of this system be 

H=^ + ^ + U{x) + V{x,y), (7.5) 
2M 2m 

where x and y are the particle coordinates (their dimensionality will be irrelevant for us, so we may consider 
them one-dimensional), U(x) is the potential energy of the heavy particle that scales proportionally to M 
as M — > oo, and V(x,y) is the potential of the particle interaction and of the light particle alone and is 
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independent of M. The wave function fy(x,y,t) of the system satisfies the Schrodinger equation 

We take the point in (x, ^-configuration space where the system is at initial time t = 0, as a coordinate 
system's origin. Then the initial action function p{x, y) is a power series in x and y, and it may be presented 
as a sum p(x, y) = Pm( x ) + Pm(x, y), where pu collects all the terms of the series with the powers of x alone, 
and p m the remaining terms, which contain nonzero powers of y. The wave function 9(x, y, t) may always be 
represented as a product A{x, t)e^' h, ' S ( x ' t ' <f>(x, y, t) with real functions A(x,t) and S(x,t). We have then from 
Eq. (HU) 



^ ,'h h 2 8 2 » . i S 



Kd__tf_j¥_ v \, (h 1 dS h 2 1 dA\ aft ft 2 d 2 </> 
ldi~2m'dy 1 + J + \~i M ~dx ~ M A~dx~ J dx ~ 2M dx 1 



(7.7) 



Let now the function A(x, t)e^^ s ^ x '^ cancel the first term in (|7.7j) . i.e., it satisfies the equation 

'h d Ti 2 d 2 \ . i , „ 

with initial condition A(x, 0)e^^' s ^ x,0 ^ = e^ h ^ PM ^ x> . The function <ft(:r, y, t) must then cancel the second term 
in (|7.7p . i.e., satisfy an equation 

To investigate the A/ — > oo limit, expand the functions 5, A, and in powers of 1/M as 

fe=0 fc=0 fc=0 

where 5^(a;,t) is proportional to A/ while coefficients (x, t), A^ k '(x,t), and (j)^ k \x,y,t) are M- independent, 
and neglect all contributions with positive powers of 1/M. For S^(x, t), we have then the classical Hamilton- 
Jacobi equation 

dS^ 1 ^dS^\" 

' +U = 0, (7.11) 



dt 2M \ dx 

and for A^°\x,t), a continuity equation 



dAi ° }2 1 9 V> 2 ^)=0. (7.12) 



at A/ dx V da; 

The solution of Eq. (|7.11[) is given by integrals of the Lagrangian function along classical trajectories in the 
potential U(x), and so S^ c \x,t) will scale proportionally to M as M — > oo, as expected. The velocity of the 
heavy particle will converge for M — > oo to an Af-independent limit v(x,t) = (1/M) dS^ c \x,t)/dx, and since 
the action S^ c \x,t) satisfies the Hamilton- Jacobi equation, this particle will exhibit a classical motion in the 
potential U{x). Let now x(t) be the trajectory of the heavy particle. Since it represents a macroscopic object, 
this trajectory is directly observable and, as such, known. The behavior of the light particle is described 
by momcntums p a with multi-indices a that include y at least once. These momentums are derivatives of 
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{Ti/i) In ^(x, y, t) taken at x = x(t). That means that the light particle is described by the wave function 
A(x,y,t) at a point x(t), i.e., in the M —> oo limit, by the function ip(y,t) — A^- ' (x(t), y, t). From (|7.9p . the 
function A^- ' satisfies the equation 

7 at " + v > ( } + T W ^*)-Sr = °" (7 ' 13) 

Combining the last term in 1 3[) with the first one, and letting W(y,t) = V[x(t), y, t) , we then obtain the 
equation for ip(y,t): 

h dib h 2 d 2 ib „ r , . , ,„ ,s 

iI-^w + WM ^ = °> (7 - 14) 

which is the Schrodinger equation for a light particle in the potential W(y, t) created by a heavy particle 
moving along the classical trajectory x(t). Thus in PAQD, the experimentally observed separation of reality 
into a macroscopic world that behaves classically and a microscopic one that exhibits quantum behavior in 
a classical macroscopic environment is not postulated as in standard quantum mechanics, but obtained as a 
direct consequence of its equations of motion. 



8. The theory of quantum measurements 

Besides different equations of motion, the difference in the measurement procedure is probably the most im- 
portant difference between classical and quantum theory. For every physical quantity, quantum mechanics 
specifies a corresponding linear hermitian operator O. In 'PAQD, Oip{x.)/ip(j(.) may be interpreted as a nu- 
merical value, which this quantity has if a particle with wave function ip happens to be at a point x. If ip 
is an eigenstate of O, then this value is the same for all x (i.e., for all possible trajectories of the particle) 
and is a corresponding real eigenvalue of O. If, on the other hand, -0 is n °t an eigenstate, then this value 
will be different for different x, and for a given point x will in general be an arbitrary complex number that 
would have been the result of a measurement of O, if this measurement had its classical meaning. In quantum 
theory, however, the situation is more complicated. Indeed, in contrast to classical theory, which deals with 
macroscopic objects, quantum theory describes microscopic ones, whose properties are usually not directly 
observable. In order to find the value of any physical quantity that such objects possess, one has to produce 
the interaction of this quantity with another one that is observable, and to infer the value of the quantity 
of interest from the reaction of that observable quantity. The observable quantity is a characteristic of the 
"apparatus" , and may have a macroscopic character, like the position of a pointer, or a microscopic one, as in 
a Stern-Gerlach experiment, where the measured quantity is a particle's spin and the observable (or rather, in 
this case, detectable) quantity is this particle's position, and the role of an apparatus is played by the particle 
itself. Thus the measurement procedure in the quantum domain is highly indirect, which causes its peculiar 
properties. To analyze them, we will apply our theory to the combination of a particle and an apparatus. We 
will identify several different kinds of quantum measurements, and consider them in turn. 

8.1. von Neumann's measurements with discrete spectrum 

The measurements of the first kind were originally investigated by von Neumann |32j . and so we will call 
them von Neumann's measurements. In this subsection we will consider the case where the spectrum of a 
measured observable O is discrete. According to von Neumann, if the apparatus performs a measurement of 
this observable, and the particle's state is its eigenstate ipi (which is assumed to be normalized, J \ipi\ dx = 1), 
corresponding to an eigenvalue Oi (so that, in "PAQD, the quantity O has the value of Oi for arbitrary 
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position x of the particle) then the reading of the apparatus should have the corresponding i-th value, clearly 
distinguishable from others. In more detail this means the following: Before the measurement, at initial time 
t = 0, the apparatus is set into the state <p(y), where y is the apparatus coordinate, which is assumed to 
be directly observable. We also assume that <p(y) is centered at y — and has width Ay. Since before the 
measurement a particle and an apparatus are independent, if a particle is in a state -0^, then an initial wave 
function of the combined particle-apparatus system is (p(y)ipi(x). If (y,x,t) is the result of an evolution of 
this state during the measurement, then it is required that for t larger than the duration of measurement At, 
^^(y,x,t) should be centered around some yi(t) and have such a width that the overlap of ^W(j/,x,t) and 
\&")(y,x, t) in y-space could be neglected for all j ^ i (in PAQD, vpW an d $W are analytic functions, and 
so they always overlap, but we can require each of them to be negligibly small in the area where the other 
one is centered). In other words, over the measurement time At different packets should diverge in y- 
space far enough to make their overlap negligible. By observing the value of y after the measurement, we 
can then infer the value of O before it. In particular, if the measurement time At is so short, and particle- 
apparatus interaction Hamiltonian Hi nt is so strong, that during the measurement all other terms in the total 
Hamiltonian of the combined particle-apparatus system may be neglected compare to Hi ntl and if Hi nt is 
proportional to O, then the wave function of the system ^^(y,x, At) immediately after the measurement will 
have the form ip^ (y)ipi(x) , i.e., the particle after the measurement will remain in an eigenstate ipi of O. But 
this is not necessary. Explicit models of such a measurement are considered in |32) and, in great detail, in |33| . 

If such an apparatus is built, then an interesting situation occurs when, before the measurement, a particle 
in not in an eigenstate of O, i.e., if its wave function is ip(x) — 53i c iV'i( x ) with more than one nonzero 
coefficient Cj. We assume that this wave function is normalized, so that Y^i \ c i\ 2 — !• By the linearity of 
Schrodinger's equation, in this case the initial wave function of the combined system (p(y)ip(x) evolves during 
the measurement into ^jCi$' ! '(!/, x, f), in direct contradiction with experiment, from which we know that 
in fact the combined system will end up in one of the states ^^'(y, x, t). To save the theory, von Neumann 
postulated, besides the unitary evolution described by the Schrodinger equation, the second law of evolution, 
which acts only during the measurements: a random, unpredictable, and unanalyzable collapse of the linear 
combination ^T. Ci*S>^ (y, x, t) into one of (y, x, t) with experimentally observed probability |c, | 2 . Nobody, 
however, was able to formulate convincingly when the unitary evolution should be replaced by the collapse (or, 
in other words, what exactly allows us to qualify an experiment as being a measurement). Similar issues arise 
in other orthodox approaches to the interpretation of quantum theory. This is the essence of the quantum 
measurement problem, which found a simple and natural resolution in the framework of DBBT [121 113) . We 
now reproduce Bohm's solution of the problem using the language of 'PAQD. 

In 'PAQD, the state of the combined system is characterized by its position in configuration space and all 
its momentums, all of which evolve according to the corresponding equations of motion. As a consequence 
of this evolution, the action function, which is just a corresponding Taylor series, evolves according to the 
quantum Hamilton- Jacobi equation, while the wave function evolves according to the Schrodinger equation as 
was described above. During this process, the combined evolution of the system's position and wave function 
is such that the system normally stays in the areas of configuration space where the wave function is not small. 
Consequently, when the packets (y, x, t) start to diverge, the apparatus position y(t) will end up in the 
area where one of them, say the fc-th, is not small, i.e., near £/&(£). Now, the momentums are derivatives of the 
system's action function, i.e., (h/i)ln^2 i Ci^/^ l '(y,x,t), at the point (y(t),x(t)) (where x(t) is the particle's 
position) in configuration space, and the further the packets move away from each other the closer are these 
derivatives to the ones of (h/i) In ^ k '(y, x, t). The measurement ends when the overlap of the packets becomes 
negligible, and with it the difference between the exact momentums and the derivatives of (h/i) In ^^(y, x, t) 
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becomes negligible also. Consequently, although the wave function is still equal to Cj'J'W (y, x, t) , the motion 
of the "physical" variables, i.e., the system's position in configuration space and momentums, will be the same 
as if the wave function was equal to "J/^ (y, x, t), in agreement with experiment. This explains the apparent 
wave function collapse. The probability of observing the fc-th result of the measurement is calculated according 
to the general rules of section 6 as an integral from | ^ - Cj\&W (y, x, t) | 2 over the area where vJ/C ) (y, x, t) is not 
small, i.e., around yk(t) in y-space and all x-space, and since all (y, x, t) are normalized and don't overlap, 
this integral, again in agreement with experiment, is equal to \ck\ 2 . 

The following features of von Neumann's measurement procedure deserve special mention: 

- Although measurement statistics are determined by the wave function of the particle alone, the result of 
every individual measurement (unless the particle was in the eigenstate of O before it) is determined by the 
full VAQD states (i.e., positions and all momentums, or positions and wave functions) of both particle and 
apparatus. 

If the paticle was not in the eigenstate of O, then the measurement's result Ok is completely unrelated to 
the value Oip(x)/ip(x) (where x is the particle's position) of observable O before the measurement. This and 
the previous note mean that unless the particle was in a corresponding eigenstate, O's observed value is not 
really measured, but rather created by the particle and apparatus jointly in the process of a measurement. 
What is measured (by the corresponding relative frequencies of a series of measurements) is a set of values of 
the squared amplitudes |q| 2 . 

- Unless the interaction Hamiltonian is proportional to O and satisfies other requirements discussed above, 
after the measurement the particle doesn't have to be in a state with a definite O value, let alone the state 
with O equal to the measured eigenvalue Ok- 

- To successfully perform a measurement, the apparatus doesn't have to be macroscopic. The only necessary 
condition is that the packets (y, x, t) with different i do not overlap after some time At (the duration of 
the measurement). In a Stern-Gerlach experiment, where the apparatus is the particle itself, the measurement 
ends and the wave function collapses not when the particle is detected after passing the magnet and we learn 
the spin measurement's result, but earlier, when the wave packets corresponding to the different spins cease 
to overlap. See, however, the next note. 

- Although we are discussing the wave function collapse, the "empty" packets $W(j/,x,t) with i ^ k do not 
disappear, but just move away from the "active" packet (y, x, t), so that their contribution to momentums 
and, therefore, their influence on the dynamics of the system vanishes. If, in their future evolution, all or some 
of the packets \£W have again overlapped with 'I'W, then the measurement would be "undone", the wave 
function would "uncollapse" , and the value of O would again become undetermined (in a sense that instead of 
being equal to \I>( fe ), the wave function would become equal to the linear combination of and overlapping 
packets \&W). This overlap, however, should happen in an x-y space of dimensionality dimx + dimy. For 
the purpose of this argument, y should include all coordinates of the apparatus and its environment that are 
connected by a chain of nonnegligible interactions. Consequently, while dimy is small (as before the particle is 
detected in a Stern-Gerlach experiment) such reversion of the measurement can, in principle, be accomplished. 
However, as soon as dimy becomes macroscopically large (as when the particle is detected or observed by any 
macroscopic, conscious or not, observer) the reversion becomes practically impossible, and its possibility may 
be neglected. 

- The set of possible final states {\I>W (y ; x , t)} of the system is predetermined by the measurement apparatus 
and does not depend on the initial wave function of the particle. Consequently, after the wave function collapses 
into one of these states, all information about the particle's initial state, and all influence of this state on the 
future history of the system is lost. On a positive side, that means that von Neumann's measurements are 
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convenient for experiment preparation. Indeed, after the observable O is measured and found equal to Ok, 
say, we know that the system is prepared in the state (y, x, t), regardless of the initial wave function of a 
particle. 



8.2. von Neumann's measurements with continuous spectrum 

An analysis, similar to that just presented, is also possible when the spectrum of a measured observable O 
is continuous. First consider the case when O is not a particle's position. Here it will be easier to use an 
explicit consideration, based on a particle-apparatus interaction Hamiltonian Hi nt proportional to O. Following 
|12[ 1131 1321 133) . choose it in the form Hi nt = g(t)Op y , where p y = —ihd/dy is the momentum conjugate to 
the apparatus position y, and the factor g(t) represents the switching of the interaction on and off. Assume it 
has an impulsive character, so that g(t) — go for < t < t and g(t) = for t < and t > r, where r is the 
duration of the measurement. Consider the limit of very small r and large go. The influence of the particle's 
and apparatus' own Hamiltonians on the evolution of the wave function during the measurement may then be 
neglected compared to the influence of Hint , so that between t = and t — r the Schrddinger equation may 
be approximated by 

ih — = H int ^ = -ihgoO — . (8.1) 
at ay 

Let ij)a{x) be eigenfunctions of O, Oip a {x) — aip a {x), normalized so that 

i>* a ,(x)i> a „{x)dx = 5{a' -a"), (8.2) 



and c(a) be the coefficients of an expansion of the initial particle's wave function ip(x,t = 0) in an integral 
over them: 

t/>(x,0) =/ c(a)ip a (x) da. (8.3) 



As before, assume the initial wave function ip(y) of the apparatus to be centered at y — 0, have width Ay, 
and be normalized, f \ip(y)\ 2 dy = 1. The total wave function of a system 9(x, y,t) may be expanded as an 
integral over ipa(x) as 

V(x,y,t)= c a (y,t)ip a (x)da. (8.4) 



It is easy to see that the expansion coefficients c a {y,t) must satisfy the equation 

dc a (y,t) dc a (y,t) 

-^ = - 9 ° a ^y— (8 - 5) 
with the initial condition c a (y, 0) = c(a)ip(y). Then the solution for c a is 

c a (y, t) = c a (y- g at, 0) = c(a)<p(y ~ goat), (8.6) 

so that at the moment t = r at the end of the measurement, the system's wave function will be 

V(x,y,r) = c(a)ip a (x)<p(y - g ra) da . (8.7) 



It is convenient to introduce a new apparatus coordinate y — y/(go r ) and new function <p(z) = (p(gorz), 
which becomes negligible when \z\ > a, where the half-width a = Ay / (2gor) . The system's final wave function 
can then be written as 

*(x,y,r)= f c(a)i/j a {x)tp{y ~ a) da . (8.8) 
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The directly observable coordinate y plays now the role of a pointer for the measurement of O. Indeed, if 
the initial wave function of a particle ip(x,0) is an eigenstate of O, say ip ao (x), then c(a) = S(a — do) and 
^(x,y,r) — ip ao (x)tp(y — do). Since tp(y — ao) vanishes for y — ao| > o, the value of y at t — r will be between 
ao — a and ao + c, so that y points to the correct value of O with precision a. It is assumed, that parameters 
Ay, go, and r may be chosen at will, and so a can be made arbitrary small. Consequently, although never 
exact, the measurement of O can be made arbitrarily precise. 

In a general situation, when tp(x, 0) is not an eigenstate of O, the scaled position y evolves during the 
measurement according to the equations of motion, and ends at t = r at some yo- Similar to the discrete 
spectrum case, because of the properties of the function tp, the evolution of the system's coordinates and 
momentums will then be the same as if the wave function at t = t instead of being ^(x, y, r), Eq. (|8.8[) . was 
equal to 

^ Vo \x,y,r)= c(a)i> a (x)tp(y-a)da. (8.9) 

Jyo—a 

Since c(a) and tpaix) are smooth (analytic) functions of a, for sufficiently small a the function ^^ Vo ^(x, y, t) 
may be approximated with arbitrary precision as 

_ rvo+c 

(x, y, t) = c(y )i>y (x)f(y) , f(y)= lp(y -a) da, (8.10) 

J y -o 

so that the function &w°'(x,y,T), to which the system appears to collapse, is an eigenstate of O with the 
eigenvalue yo- The probability p(ag, day) to find the value of O (i.e., the value of yo) between ao and ao + dao 
is equal to p(ao)dao, where p(ao) is the corresponding probability density. By the general rules of section 6 we 
have for it 

p{ao) = goT J \^{x,goTa ,T)\ 2 dx = g r J \c(a)\ 2 \ip(a - a)\ 2 da , (8.11) 

where we used the normalization condition for the functions ip a {x) and the factor g$T appears because dy = 
g$r dy. For sufficiently small a, \c(a)\ 2 in the integrand may be again approximated by |c(ao)| 2 with negligible 
error, and using normalization of the function tp, we obtain, in agreement with von Neumann's postulate and 
experiment, the standard result p(ao) — c(ao)| 2 - 

To summarize, von Neumann's measurement procedure of an observable with continuous spectrum that is 
not a particle's position is similar to the one with discrete spectrum, and has the same, listed above, properties. 
In particular, unless the particle was initially in the eigenstate of O, the measurement's result is unrelated to 
the value Oip(yL)/ip(x.) of O before the measurement, and with arbitrary precision the set of possible final states 
of the particle-apparatus system is predetermined and does not depend on the initial state of the particle. For 
the measurement procedure, considered above, the corresponding set of possible particle final states is just the 
set {tp a (x)} of eigenstates of O. 

We now consider von Neumann's measurement of a particle's position, and show that, in contrast to 
other physical quantities, this measurement results in the true VAQD particle position. Indeed, we have, 
obviously, for the eigenfunctions ip a {x) of an operator O = x and coefficients c(a) of expansion of the particle's 
wave function tp(x, 0) in this case, ipa{x) — 8{x — a) and c(a) = ip{a, 0). Although 8(x — a) is not an analytic 
function of x, and so cannot be considered a legitimate "PAQD wave function, we still can use it in intermediate 
mathematical transformations. The integral over a in Eq. (|8.8[) can then be immediately calculated to give 
^(x,y,r) — ip(x, 0)(f(y — x). Let now the particle's position before the measurement be xq. Assume that the 
function tp(y) is real. Then the motion of the particle during the measurement may be neglected. The final 
value yo of y will now be such that ^{xo, goTyo, t) does not vanish. Due to the properties of the function tp, 
that means that yo should be in a er-vicinity of xq, \yo — Xo\ < c, as was asserted. The collapsed wave function 
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t), Eq. (18.91) . is easily calculated to be equal to ip(x,0)ip(y — x) for |a; — 7/0 1 < c and equal to zero for 
\x — yo| > a. It is, therefore, contained in a cr-vicinity of ya, and for a —> 0, as in other cases of von Neumann's 
measurements, loses the memory of the particle's initial state. 

8.3. Position measurements of the second kind and the double-slit experiment 

In practice, von Neumann's measurement is never used for a particle's position. The real position measurement 
is carried out by such devices as a bubble chamber or photographic plate. As we will see, wave function collapse 
and some other important features of this measurement are significantly different from those of von Neumann's, 
which justifies calling it the measurement of the second kind. 

The measurement of a particle's position by a photographic plate or in a bubble chamber may be described 
as follows. The physical state is filled with microscopic detectors (molecules of photo-emulsion for photographic 
plate, or of overheated liquid for bubble chamber), which change their state (chemical changes in emulsion, 
ionization in a liquid) if the measured particle passes in close vicinity to them. Due to the special physics of 
detectors, this changed microstate evolves then in such a way as to produce directly observable macroscopic 
changes (dark spot on developed plate, bubble in a chamber). The detectors with such changed state mark 
the position of a particle. 

In our analysis of this procedure, we again use the fact that the wave function in VAQD has its own 
dynamics, the same as in standard quantum mechanics, and so its evolution may be analyzed without reference 
to a particular particle's trajectory that is responsible for this wave function's existence, and which may be 
included in the analysis later. Then the following crude model may be suggested for the description of a 
position measurement. Consider first just one detector, fixed at a point with position x. Let the detector 
itself be characterized by the parameter y, which in the initial state is close to zero, so that the initial wave 
function of the detector is, for example, cxp(— y 2 /4) (we will not worry about wave function normalization 
here). Let the particle's wave function be ip(x), so that an initial wave function of the particle-detector system 
is ip(x) exp(— y 2 /4). Let the physics of the detector and its interaction with the particle be such that within 
a short measurement time Ai, y moves from the vicinity of zero to the vicinity of some F > 1, so that the 
detector's wave function becomes, for example, exp \—(y — K) 2 /4] , if during this time the particle's distance 
from the detector was less than some characteristic distance a. The evolution of a particle's wave function 
during the measurement time due to its own Hamiltonian, i.e., without interaction with the detector, will be of 
no importance for us. We can, therefore, consider the measurement to be instantaneous, i.e., At to be so small 
that the change of the particle's wave function during the measurement due to its own dynamics is negligible. 
The wave function of the particle-detector system immediately after the measurement can then be written in 
a general form as 



W = ip(x) < a(x — J, y) exp 



\iv-Y? 



- hi,- - 1/ ) (>\p j — i y 2 ) } . (,Y!2) 



The functions a and b reflect the physics of the particle-detector interaction. All we know about them is that 
a(x — x, y) vanishes and b(x — x, y) converges to 1 when the distance from x to x becomes larger than a, and 
that b{x — x, y) vanishes when this distance is smaller than a . The dependence on y is included in a and b for 
generality, and is supposed to leave the general character of y-dependence of the corresponding terms intact, 
i.e., the probability density is concentrated near y = Y in the first term, and near y = in the second. Now 
if there are many detectors like that, then the initial wave function will be tp(x) Y\ i exp(— y 2 /4), and the wave 
function after the measurement will be 



* = ip(x)Y[^a(x-Xuyi)exp -^(y l -Y) 2 + b(x - x it y t ) exp yfj j 



.13) 
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where j/i is the y-coordinate of the z-th detector, and X; is its position. To avoid unnecessary complications, we 
will make the simplifying assumption that a- vicinities of different detectors do not overlap and, at the same 
time, do not leave any places in the x-space uncovered. The after-measurement wave function may then be 
rewritten as 



Let now the particle's position x p at this moment happen to be in a a- vicinity of the fc-th detector, x p ~ x^. 
Then obviously in the sum over j in (|8.14p . all terms except the fc-th, will give a negligible contribution to 
momentums (i.e., derivatives of (h/i) In VP at that x p and proper j/^'s) and so the future motion of the particle 
and detectors will proceed as if the wave function was equal to this fc-th term, i.e., underwent a collapse 



In this state, the y-coordinate of the fc-th detector will then be found near Y, and all others will be near zero, and 
the particle's probability density, although far from being a delta- function centered at x p , will be concentrated 
in x p s and the fc-th detector's cr-vicinity, where a may be considered as a measurement precision. With this 
precision, therefore, a position measurement of the second kind, like its von Neumann's counterpart, measures 
the true VAQD position of a particle. We note, however, that although for both kinds of measurements the 
final selection of a member of a linear superposition, to which the wave function would collapse, is made 
by some variable which may have one, and only one value, for von Neumann's measurement this variable is 
the apparatus position y, while for a position measurement of the second kind, it is the measured particle's 
position x p . Also, for a position measurement of the second kind, the precision a is fixed by the physics of 
detectors and so cannot be made arbitrary small. Consequently, while for an arbitrarily precise von Neumann 
measurement, the final wave function becomes equal to one possible function from the predetermined set of 
them, the final wave function after a position measurement of the second kind does depend on the initial wave 
function of the particle. Indeed, its x-dependence is essentially given by the product ip(x) a(x — xt, Y), i.e., is 
equal to the initial function ip(x) modulated by a factor a(x — Xk, Y). If this factor is smooth enough, and if 
the characteristic wavelengths in ip(x) are much smaller than <r, then the packet ^){x)a(x — Xk,Y) will keep 
propagating without spreading much along a trajectory that is close to the one the particle would have by 
itself, i.e., if its position was not measured. In a bubble chamber, this packet will then trigger other detectors, 
thus producing a track which approximates the particle's unperturbed trajectory. 

The same consideration may be also applied to the double-slit experiment discussed by Feynman [15] . In 
this case, variable x in Eqs. (|8 . 1 3[) - (|8 . 1 5[) will denote the coordinate on the screen, and Xi, i = 1,2 — the 
position of the i-th slit. Without detectors, the wave function of the particle immediately behind the screen 
would be equal to ip( x ) Hj=i i a ( x ~ x j)> wnere ^{x) is the wave function in front of the screen, and the 
"shadow function" a(x — Xj) is nonzero only for x inside the j-th slit. Propagating away from the screen, the 
waves from the two slits would overlap and create an interference pattern. On the other hand, in the presence 
of detectors the wave function will be given by Eq. (|8.14j) with indices i,j there taking the values of 1 and 2. 
The condition that detectors work well and allow to determine through which slit the particle have passed 
means then exactly that the packets from the two slits remain well separated with respect to coordinates y\ 
and ?/2, and so the interference between them is impossible. As was explained above, if the particle have passed 
the slit fc and was detected there, then its future motion and the motion of detectors will be the same as in the 
state with the wave function (|8.15p . i.e., as if the other slit was closed. In agreement with [TS], the observation 
of an interference is, therefore, incompatible with the detection of the path chosen by the particle. These two 




(8.14) 
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operations are just mutually exclusive: the interference happens when the packets overlap, while the detection 
of the path requires them to be well separated. Note that this conclusion remains perfectly valid even when 
the detectors are microscopic, like the one-bit detectors discussed in (54"] . 

9. Nonlocality, analyticity, and covariance 

Although in PAQD, as in classical mechanics, particles move along well defined trajectories, the equations of 
motion in these theories are fundamentally different. The only momentums that contribute to the equations of 
classical mechanics arc the first order momentums pj k , where k runs from 1 to the number of particles n p , and 
for every k, jk runs from 1 to the dimension of physical space n. Every momentum pj k , therefore, is "bound" to 
a corresponding particle k, and changes only due to the presence of forces, described by a potential function U . 
In nonrelativistic mechanics, the forces corresponding to this potential normally vanish with distance, while 
in relativistic cases the potential propagates with finite speed, which is not larger than the speed of light c. 
Classical mechanics is, therefore, local: to predict the behavior of a particle separated by a large distance from 
others during some time At, one doesn't have to know what happens further than the distance of cAt from 
it. This locality, we see, is a consequence of the fact that in classical mechanics particles influence each other 
only through the action of the potential, which has the described properties. 

In PAQD the situation is different. To avoid tedious manipulations with a multiparticle Hamiltonian (I7.2[) . 
we may simply make all masses equal to each other and denote them as m. The Hamiltonian (|7.2p will 
then look exactly like the one-particle Hamiltonian (|3.5[) . but in an n p n, rather than in an n-dimensional 
space. Correspondingly, Eq. (|3.24[) for the evolution of momentums will hold, with summation over repeating 
indices j there running from 1 to n p n. Now if particles are entangled, i.e., the system's action function is not 
equal to the sum of separate particles' actions (or system's wave function to the product of separate particles' 
wave functions) then there exist nonzero momentums p^ with "mixed" multi-indices /i, which include indices 
from different particles. Eq. (|3.24p will then interconnect the time evolutions of all possible momentums p M , 
and with them of particles' velocities. Since the momentum-dependent part of (|3.24p does not depend on the 
particle positions, and all momentums are taken at the same time, they obtain the status of global variables: 
each momentum affects the time evolution of all others at the same moment of time, independently of the 
particle positions and the distances between them. Thus in this new (i.e., nonclassical, "nonpotential" ) way, 
the particles in PAQD influence each other on the whole hypersurface t = const instantaneously, and over 
arbitrary distance. Clearly, the reason for this nonlocality is that an analytic function is a fundamentally 
nonlocal object — the set of its derivatives in any point of space determines its behavior arbitrarily far from 
this point. 

The nonlocal kind of behavior described above, is, according to Bell's theorem, necessary for any theory 
that dynamically derives experimentally observed nonlocal correlations between entangled particles, rather 
than just predicts them, as does standard quantum mechanics [351 [35]. "PAQD is built as an "ODE side" of 
quantum mechanics, which always agrees, of course, with its "PDE side", i.e., the Schrodinger equation and 
the conventional theory based on it. As such, "PAQD must be nonlocal: if it were local, so would the standard 
quantum mechanics. Note also, that "PAQD does not conflict with our intuition: indeed, our intuition is 
classical, but the classical limit of "PAQD is just the usual, completely local classical mechanics! PAQD thus 
has the desired feature of being a fundamentally nonlocal theory with a local classical limit. 

The described nonlocal behavior was first discovered in the framework of DBBT and discussed extensively 
there [13] , It was soon realized that DBBT's nonlocality is in perfect accord with the requirements of Bell's 
theorem and is, in this respect, welcomed |35) . There remained, however, a difficult question about the 
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theory's relativistic invariance. The influence, propagating with infinite speed, seems to be in an obvious 
conflict with the requirements of special relativity theory. This concern is addressed in |13j , where it is proved 
that such influence cannot be used for transmission of superluminal signals. Still, there is the other concern: 
propagation of influence with infinite speed requires a selection of preferred reference frame, in which this 
propagation happens along surfaces t = const, in contradiction with the spirit of the theory of relativity, which 
demands that physical laws must be the same in every inertial frame of reference. This is generally considered 
to be a serious problem for DBBT [37]. We will now show, that, thanks to the additional requirement of 
analyticity, "PAQD may be formulated in an arbitrary analytic foliation of space-time, and will have the same 
form in each of them. Our consideration will be nonrelativistic. It will be argued at the end of this section, 
however, that its relativistic version, although still nonlocal, will be not only Lorentz invariant, but can be 
also made generally covariant. 

Indeed, consider an arbitrary analytic foliation of space-time, generated by a single-valued analytic func- 
tion /(x, r), i.e., a partition of space-time into 3-dimensional hypersurfaces (leaves of foliation) 

i = /(x,r), (9.1) 

where r parameterizes hypersurfaces (we can, for example, conveniently require r — /(0,r)) monotonically, 
so that 9/(x, t)/8t > for all x and r, and such that the whole space-time is covered (so that every point 
(xo,£o) belongs to some hypersurface, i.e., to = /(xo,to) with some To). In relativistic theory we require the 
surfaces r = const to be space-like. We will call this foliation /-foliation. The standard partition of space-time 
into surfaces t — const ("standard foliation") corresponds to a function /(x, r) = r. We now want to introduce 
wave functions, defined on surfaces r = const, rather than t = const. For a one-particle case, the wave function 
ip(x,t) was introduced as a solution of a Schrodinger equation, analytic with respect to x for every t. It is 
then also analytic with respect to both x and t, and so the function ^^(x, r) = ipix., /(x, t)) is analytic with 
respect to x and r. To define a wave function on surfaces r = const in a multi-particle case, we will borrow 
from relativistic theory the multi-time formalism [38] . where each particle has its own individual time, and 
the multi-time wave function of n p particles ipfei, t\, . . . ,x„ p ,t nj) ) depends on positions and times of all of 
them. Detailed analysis of the physical meaning of this wave function and of corresponding analytical quantum 
dynamics will be a subject of relativistic consideration. In nonrelativistic theory, where the interaction between 
particles is mediated by an instantaneous potential function, the multi-time formalism can be defined only for 
particles that do not interact with each other (but can interact with an external potential) . It will be sufficient 
for our purpose, however, to consider such noninteracting particles, because here we are only interested in 
nonlocal correlations, caused by entanglement, and not in correlations due to an interparticle interaction. The 
multi-time wave function then satisfies the system of equations 

ih^ H k ) il>(xi,ti,...,x np ,t np ) = 0, k=l,...,n p , (9.2) 



dt k 



where 



h 2 

H k = -- A k + U k (x k ,t k ) (9.3) 

2m k 

is the Hamilton operator of the fc-th particle in the external potential U k (n,t), A k being a Laplace operator, 
acting on the coordinates of the fc-th particle X&. The wave function ip{x.i,tt, ■ ■ ■ ,x„ p ,£„ p ) may be obtained 
by path integration over all the paths such that for every fc = 1, ...,n p , the paths for the fc-th particle 
terminate in a point (x k ,t k ). The action functions S(x.t,ti, . . . ,Xn , t n ) and i2(xi,ti, . . . ,x n , t np ), defined 
from iplpci, t\, . . . , x„ p , t nv ) as in Eq. (|7.I I) , satisfy a system of quantum Hamilton- Jacobi equations 

— +H% = Q, — +H * = 0, k = l,...,n p , (9.4) 
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similar to Eq. p.4[) , with H% and Hj} given by Eqs. (|3.6I) and (|3.7[) . where the index j in those equa- 
tions means the derivative with respect to the j-th coordinate of the fc-th particle, and the potential U is 
understood as Uk(xk,tk)- As will become clear soon, it is appropriate to postulate, in a straightforward gen- 
eralization of Eq. p. 81) and corresponding one-time theory, that if in a state described by a wave function 
tp(xi,ti, . . . , x n , t np ), the particles have space-time positions (xi,ti), . . . , (x n ,i n ), then their velocities are 
given by 

dH s dH R 1 

v fe = ^ = ^ = — S Xfe (x l5 x np , i n J , (9.5) 

where S Xfc (xi , *i , . . . ,x„ p ,i„ p ) = <9S(xi,ii, . . . , x Ilp , t np )/cht k , and similarly for i? Xfc . 

With so-defined multi-time wave function, the wave function on a hypersurface r = const of any /-foliation 
is obtained by placing each particle on this hypersurface, 

■ • ■ ,x„ p ,r) = ^(xi,/(xi,r),...,x„ p ,/(x„ p ,r)) , (9.6) 

and it is an analytic function of all x^ and r. The whole theory developed above for a standard folia- 
tion may then be reproduced for an arbitrary analytic /-foliation. The T-evolution of the wave function 
ipW (xi , . . . , x„ p , r) is governed by the equation 

ih ^!l = H(f)^fl (9.7) 

OT 

The transformation of space coordinates does not affect our nonrelativistic analysis, and so we will use the 
same coordinates in all foliations. Then 

n p 

HU) = g £(/) , njp = Mx k ,r)H k , (9.8) 

k=l 

where f T denotes the derivative of / with respect to its second argument, i.e., f T (xk,r) — <9/(xfc, t)/8t. In 

what follows, we also need the derivative / x of / with respect to its first argument, / x (xfe, r) = df(x.k, T)/dxk- 

The action functions p^ (xi , . . . , x„ p ,r), S W (xx , . . . , x„ p , r ) , and B.W (xi , . . . , x„ p , r ) on /-foliation are again 

obtained from ipW as in Eq. (|7.1[) . and then the quantum Hamilton- Jacobi equation, the momentums p^J\ 
( f) ( f) 

Sa , and RJ and the equations of motion for them are introduced in the same way as for a standard foliation. 
In particular, the quantum Hamilton- Jacobi equation for S™> is 

^+i^) S =0, (9.9) 

OT 

where 

rip 

H Uf = E H k )S > H i ff = M*k,T)Hg , (9.10) 
k=i 

and particle velocities in the /-foliation, i.e., with respect to a new "time" r, are given by the usual relation 
v[ f) = dH^ S /dSi{\ where S { JJ = dS^/dx k , or, using (l9~T0l) . 



,(/) 



f)H s 

frto'T)-^. (9.11) 



Note, that the derivative over x^ in Sx k is taken along the leaf of the /-foliation, i.e., for r = const, contrary 
to the derivative in S Xk , which is taken for t — const. 
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We can now show that velocities v!l , Eq. (|9.11[) . and v^, Eq. (|9.5[) . correspond to the same motion of the 
fc-th particle. We note first that if this particle moves from point Xfc on leaf r to point X& + dxj- on leaf r + dr, 
so that its r- velocity is = d^k/dr, then by Eq. (|9 . 1 1) we have for a corresponding time interval 

dt k = f T (xfc , t) dr + f x (x fc , r ) dx fc 

Consequently, t- velocity v fc = dx k /dt k should be equal to v£ /(/ T + /x v i ) 

v«=v fc (/ T + / x v^). 



or 



(9.12) 



(9.13) 



To demonstrate that this relation between and v& does indeed take place, we need to express S Xk in 
Eq. (JUSJ) through SiV in Eq. (|9H|) . For the space-time of the fc-th particle, consider the surface r = const, 
or tk = /(xfc, t). We have for the derivatives of the action function £(xi,ti, . . . , x„ , t n ) along this surface 



as 



But 95/9x fe | T=const = SL{\ and so 5 Xfc = S^J + / x i?f , and therefore 



as 



Xfc —const 



9x fe 



(9.14) 



= S Xfc -tff/ x (x fc ,r) 



(9.15) 



Now multiplying this equation by / T (xfc,r) and using Eqs. (|9.5j) . (|9.1ip . we immediately obtain the desired 
relation ([9TT3]l . 

Thus, for any function /, the description provided by a foliation-based wave function ij}^> agrees with the 
one provided by a multi-time wave function ip. We have, therefore, the following situation. In every foliation, 
the state of a system is described by particle positions and by all momentums. Being the derivatives of the 
logarithm of a wave function along the leaves of the foliation, momentums depend on the foliation chosen, 
and so in any given system's state, there are different sets of momentums, corresponding to different possible 
foliations. As for a standard foliation, by the equations of motion momentums, corresponding to any foliation, 
are global variables — they are bound to the leaves of their foliation, rather than to the points on these leaves, 
and influence each other (and, consequently, the particle velocities) over the whole leaves of this foliation 
instantaneously. 'PAQD, therefore, does not require a preferred frame of reference: in every foliation, the 
theory, expressed through the foliation's momentums and Hamiltonian function, looks the same. At the same 
time, considerations, based on different foliations, agree with each other in terms of actual motion of particles, 
because they all predict the same motion as the consideration based on a common object, the multi-time wave 
function ip(pc±,ti, . . . , x n , t n ), as was just discussed. On the other hand, this multi-time wave function in all 
_R 4 " p may be uniquely obtained, for any / and To, from a function ipr{ (xi, . . . ,x„ ) = t/>(^)(xi, . . . ,x„ ,to) 
by path integration or by solving equations (|9.2[) . and so every function i\r TQ contains the same information as 
the multi-time wave function ip in all R 4rlp . 

In relativistic theory, the leaves of foliations corresponding to different Lorentz frames are flat, and the 
angles between different foliations' leaves correspond to relative velocities of respective frames. In addition to 
a time transformation, a Lorentz transformation of the space coordinates inside the leaves should be done. 
It seems then reasonable to expect that in relativistic theory the set of momentums, corresponding to each 
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Lorentz frame, will behave as described above, i.e., momentums will affect each other instantaneously in this 
frame, and in each frame the theory will be the same. Also, if there are two space-like separated entangled 
measurements, then neither of them can be considered as causing the result of the other. Indeed, in different 
foliations their time order will be different, and their results are unambiguously determined by the state 
(i.e., positions and all momentums, or positions and wave function) on any leave of an arbitrary foliation. 
Foliations with nonflat leaves (i.e., leaves that in different space-time points are angled differently with respect 
to the leaves of "inertial" foliations) will appear, when the theory is formulated in reference frames with local 
accelerations. Additional terms corresponding to inertial forces, proportional to particle masses, will then 
appear in Hamiltonian functions, so that every foliation will come with its own field of these forces. The 
equivalence of different foliations, i.e., the general covariance of the theory, can then be restored in a standard 
way by introducing a gauge field, which would adsorb the potential of inertial forces, in what seems to be a 
natural route leading to a gauge theory of gravitation [39 . 

10. Analytical quantum dynamics of particles with spin 

In this section we show how to describe in "PAQD particles with spin. Since the wave function of a particle 
with spin s is a 2s + 1-component spinor0 there seem to be two possible ways to include spin in the theory. 
The first one is to somehow define corresponding 2s + 1 complex or 2(2s + 1) real action functions, satisfying 
evolutionary equations of (|2.1I) type with Hamiltonians that depend only on derivatives of the action and satisfy 
Hamiltonian conditions. This, however, docs not seem to be possible. Indeed, for the Schrodinger equation 
we passed from the wave function to its logarithm, the action function, in order to obtain an evolutionary 
(namely, quantum Hamilton- Jacobi) equation with Hamiltonian that depends only on derivatives of the action 
function, rather than on this function itself. For a multi-component wave function, this simple trick will work 
only in a trivial case when every component satisfies its own equation, independent of others. Moreover, in 
case of several, say n p , particles, one would have to find not 2s + 1 complex action functions, which satisfy 
equations of the required form, but (2s + l)" p of them! Clearly, this approach doesn't appear promising. The 
second possible approach is to transform a system of 2s + 1 equations for spinor components into an equivalent 
equation of the desired form for one new wave function. This can be done by using spin coherent states, and 
this is the approach that we will employ here. 

Spin coherent states are defined with the help of a spin rotating operator that rotates the spin state through 
an angle a about direction n. The explicit form of this operator is exp (— ians), where s is a spin operator in 
units of Ti. The rotation TZ(x, 6, ip), corresponding to Euler angles Xi &i ft is obtained as a rotation through 
the angle x about axis Oz, followed by rotation through angle 9 about axis Oy, followed by another rotation 
through angle cp about Oz, and is described by the product of the three corresponding operators: lZ(x, 0, ip) = 
e -ivsz e -i6s y e -ixs z _ \s,m) be a standard eigenstate of the spin operators: s 2 \s,m) — s(s + l)|s,m), 

s z \s, m) — m\s, m). Then the spin coherent state is defined as the maximally polarized state |s, s) rotated by 
the operator TZ(x, 0, f)'- \x> ®i = e~ l{ps * e~ lSsy e~ lxs * \s, s). The explicit representation of the coherent state is 



s 



s+m 1 ,s—m 




(10.1) 



where complex parameters u and v are defined as 




(10.2) 



7 We will denote particle's spin by the small letter s, to distinguish it from the real part S of the action function p. 
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Clearly, if u = u\ + iu2 and v = v\ + W2, where 1*1,2 and Wi,a are real, then u| 2 + \ v\ 2 — u\ + u\ + +v\ = 1, so 
that parameters u and v (or u\, U2, V\, V2) live on the three-dimensional unit sphere S 3 in the four-dimensional 
real space i? 4 . Denote the set (x, 9, <p), or corresponding sets (u,v) or (ui, U2, wi, W2), as f2, and the coherent 
state (jlO.ip as The expression (|10.1[) for it may be easily derived, for example, by using the Schwinger 
bosons representation of spin operators [3D] . It is well known (30J IH] that the system of spin coherent 
states is not orthogonal, overcomplete, and allows a resolution of unity in a Hilbert space of states with spin s: 



2s + 1 



dfl\n}(n\ = ^2\s,m}(s,m\, (10.3) 



where 

dSl = - sm6dxd6dip (10.4) 
8 

is the area element on S 3 . If wi , W2 , W3 , W4 are Cartesian coordinates in R 4 that are related to angular 
coordinates \, 9, (p and radial distance r by 

Wl,2 = TOl,2 , ti>3,4 = rvi t2 , (10.5) 

then the integration measure in R is related to dfl by 

dw\dw2dw'idw4 = r 3 drd£l, (10.6) 

which follows from the expression (most easily obtained by direct calculation with Mathematical for the 
jacobian of the transformation (|10.5I) 

d(w 1 ,W2,w 3 ,w 4 ) 1 3 . 
det — — — r — = -r sm9. (10.7) 

The action of the spin operators s± = s x ± is v and s z on spin coherent states is described by equations [40] 

s+ \n) = «fl u |n>, 

s z |^) = - (u9„ - w9„)|fi). 

Using spin coherent states, the one-component wave function, corresponding to a spin state l^) with spin s, 
is defined as a scalar product "0(^) = (^IVO- It is clear from the resolution of unity (|10.3[) that using its wave 
function, the state may be expanded over spin coherent states as 

|V-') = J d(lip(Sl)\n) , (10.9) 

so that all information about the state is contained in its wave function and vice versa. The argument O in 
the wave function describes a rotation with respect to Cartesian coordinates in three-dimensional space, and 
so when the coordinate system itself is rotated the wave function transforms accordingly. The rotations are 
elements of the three-dimensional rotation group SO (3), where the spin wave function is defined. As is well 
known, SO(3) is not simply connected: its fundamental group is cyclic group of order 2. As was discussed 
in section 2.5 and section 3, this means that the spin wave function defined in SO (3) may be double- valued, 
as it indeed is when spin s is half-integer. The universal cover of SO(3) is the group SU(2) that covers 
SO(3) two-to-one, and so the spin wave function is single-valued in SU(2). Elements of SO(3) and SU(2) 
are parameterized by points on the sphere S 3 considered above [3T], and the area element <if2, Eq. (|10.4[) , is 
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the Haar measure of these groups, so the spin wave function may be considered as defined on S 3 . Now if 
\ip) = J~) m —_ g ^m( x i i)l s i m ) is a state of a particle with spin s, then the corresponding wave function is 

S — S~\-TTh — S — 771 

tf(x, n, t) = ^s)\ J2 t, = w *) • ( 1(uo ) 

„^ s V(s + m)!(s-m)! 

Thus a wave function of a particle with spin s is defined in configuration space Q = R 3 x S 3 and is an analytic 
function of u and v, selected from arbitrary analytic functions of these variables by the condition of being a 
homogeneous function of power 2s, i.e., by condition 

(ud a + vdy)i>(x, O, t) = 2s^(x, O, t) . (10.11) 

In the space of analytic functions of u and v, an operator (l/2)(u<9 s + w9 s ) plays, therefore, the role of a total 
spin operator. From Eq. (|10.8[) we obtain the action of spin operators on a wave function: 

(fl\s + \ijj) = udv*l>(x,tt,t), 

<O|s_)-0) = vd u tp(x,D,,t), (10.12) 
(Cl\s z \t/j) = -(uda-vdv)il)(x,il,t), 

where we used that operators s + and s_ are hermitian conjugates of each other. Note, that spin operators make 
a complete set of non-trivial first order differential operators that leave a wave function in the form (|10.10|) . 
The only remaining operator, udu + vdy, gives 2s by Eq. (|10.11[) . 

The behavior of a particle of charge e in an electric field with scalar potential Aq and magnetic field B 
with vector potential A is described by 



1 e 2 



V - - A + eA - 7 Bs 

2m V i c 



(10.13) 



where c is the speed of light. For a particle with spin 1/2, Dirac theory gives for a constant 7 the value 
of eh/mc. If ip is the just-defined wave function in a spin coherent state representation, then the term Bs 
expands as 

Bs = i (B+S- + B_s+) + B z s z 

= - [B+vdu + B-udy + B z (udu - vd<o)] , 
where B± — B x ± iB y . Using matrix notations and standard Pauli matrices, we have then 



(10.14) 



„ 1 / B z B- \ f d u 

= \{u,v)K<r( d d * 



(10.15) 



As in the spinless case, expressing the wave function as 

#cn,t) = e*p(irtx,n,t)) > p(x, q, t) = s(x, n, t) + j r(x, q, t) (10.16) 

introduce the action function p(x, Q, t) and its real and imaginary parts S^x, SI, £) and —hR(x,Q,t). It is 
convenient to use a gauge divA = 0. Momentums that correspond to spin variables, such as u, u±, or \ (i- e -) 
partial derivatives of the action with respect to these variables) will be denoted by the corresponding indices. 
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The derivatives with respect to complex variables are understood as in Eq. (|2.58p . Then substituting (|10.16[) 
into (|10.13[) . obtain for a particle with spin a quantum Hamilton- Jacobi equation (|3.3p with Hamiltonian 
function 

H =L (» - l A -) 2 + eA ° + L*» - 1 < i ' f '» B " {";,)■ (iai7) 

The Hamiltonian (|10.1T|) is of the first order with respect to the spin variables. Consequently, HC1 is satisfied 
for it automatically, while HC2 is satisfied due to analyticity, as was discussed in section 2.6. The whole theory 
of section 2, therefore, is applicable, but this time in configuration space Q = R 3 x S 3 , so that at any time 
the particle has its space position in i? 3 , "internal" SU(2) position on S 3 , and all corresponding momentums. 
To guarantee that the particle's spin is equal to s, S 3 positions and momentums should satisfy 

\ (upu + wp s ) = 2s , (10.18) 
n 

which follows from Eq. (|10.11[) . Since spin operators change only the projections of spin, and not its value, 
for any Hamiltonian that, as in Eq. (|10.13|) . depends only on spin operators, it is sufficient if this condition 
is satisfied at the initial moment of time. According to the general theory of section 2, particle velocity in 
physical space R 3 is given by Eq. (|2.66[) . i.e., 

^ = ^-(p j +P j )-—A j , (10.19) 
2m mc 

while SU(2) variables evolve (see section 2.6) according to 

( *-Hls'fD = -§ (a ' i ' )B<T <io ' 2o) 



£W " |. (10.21) 



or, after hcrmitian conjugation, 

(., _ 

v J 2h \ v 

This is an equation of spinor rotation with angular velocity u> = —(■y/h)H. The time evolution of spinor (u, v) T 
(where T indicates transposition), composed of SU(2) coordinates of a particle, is, therefore, very simple: at 
any moment it rotates with this angular velocity, B being the magnetic field at the current particle's position. 
Using Eqs. (|10.20|) . (|10.2ip . it is easy to demonstrate that the value \u\ 2 + \v\ 2 is conserved along the spinor's 
trajectory, and so remains equal to one, if it was equal to it initially. 

Equations of motion for momentums couple all kinds of them: "space" momentums, with multi-indices 
composed of x, y, and z, "spin" momentums with multi-indices composed of u and v, and "mixed" momentums, 
with multi-indices composed of both kinds of variables. These equations decouple if the magnetic field B is 
spatially uniform, and the initial wave function factorizes in the form ?/)(x, f2, 0) = ijj^ x '(x, 0)^^ (0, 0) or 
p(x, Q, 0) = p( x '(x, 0) + (O,0), where the indices x and s mark the space and spin parts. The mixed 
momentums then remain equal to zero and the wave function remains factorized at all times. The space part 
of the wave/action function satisfies the equations for a spinless particle, and so the particle moves in the 
physical space as if it didn't have any spin. Using Eq. (110. 20j) . the equation 

for a spin part of an action function may be written in the form p^ — 0, where p^ = dp^ /dt + ud a p^ + 
vdvp( s \ As for every homogeneous PDE of the first order the solution p( s \ therefore, remains constant 
along the equation's characteristic curve, i.e., along the trajectory (u(t),v(t)) in S 3 described by Eq. (|10.21l) . 
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Along with p( s \ the wave function ^( s ) also remains constant, i.e., (u(t), v(t)\ip^ (t)) = const. Consequently, 
like a spinor (u,v) T , the spin part ip^ of the wave function rotates with angular velocity uj = —('j/h)G, 
exhibiting the well-known spin precession in a spatially uniform magnetic held. 

Like a spinless Schrodinger equation, Eq. (|10.13[) may be obtained from a variational principle 

3 

C dnC«.«) Y[ dx j dt = (10.23) 
with Lagrangian density 

C = Y% ^* dt ^ ~ ^ dtr) + 2^ { d ^* + Yc A ^*) i 9 ^ ~ ¥c A ^) + eA ° r ^ ~ ^* Bs ^ ■ ( 10 - 24 ) 

The spin part of the integration measure in Eq. (|10.23j) is dCl^ u ' v ^ = du\ A du 2 A dui A dv2 — — (l/4)d« A du A 
dv A du, and integration over u\, U2, v%, V2 runs over the whole space i? 4 , so that this variational principle 
defines Eq. (|10.13[) in the whole space i? 4 , and not only on the unit sphere S 3 . The phase invariance of the 
Lagrangian density C leads, by Noether's theorem, to the corresponding conservation law, which now has the 
form 

d j° + clivj + du.f + d,f = , (10.25) 

where the components of the current are given by Eq. (|5T2]l . but with i = 0, . . . , 3, u, v this time. Substituting 
there A = (i/K)tp, A* = -(i/h)ip*, A 1 = 0, and using Eq. (|10.15[) for Bs, one gets for the current 



(/,].. /"•./*') = (JV>| 2 , 2^W>*V^-^n-^AH 2 , -g|^| 2 (S,«)B<T 
= \ip\ 2 (l,v,u,v). 



(10.26) 



Note, that since this current doesn't have a radial component in i? 4 , the conservation law (|10.25j) is satisfied 
on every sphere with the center in the origin there, including a unit sphere S 3 , where we need it. The results of 
the previous sections can now be immediately generalized to the case of particles with spin. Most importantly, 
\ip\ 2 becomes the probability density in configuration space R 3 x S 3 with respect to a measure gK^I1^=i dxj, 
where dfl is a measure () 1 .4(1 on S 3 , and the measurement of spin-related physical quantities is described by 
the same theory of section 8 as for space-related quantities. 

Although, as was discussed in section 2.6, the above derivation in complex coordinates u, v is equivalent 
to the one that uses coordinates u\,U2,v\, V2, it may be instructive to present a direct derivation in these real 
coordinates. For that, it is convenient to present the last term in Eq. (|10.13[) in the form 

- 7 BsV> = - (U*du + V*dM , (10.27) 
1 ' 

where we introduced U* — U\ — 1U2 and V* = V\ — 1V2 for which, by Eq. (|10.14p . we have 



u * = In 





-V 


+ B z u)] = 




- ByVl 


+ B z u 2 ) , 


3f[i(B H 


-V 


+ B z u)} = 




+ ByV 2 


+ B zUl ), 


R[i(B. 


-U 


-B z v)] = 




+ ByUl 


- B z v 2 ) , 


$S[i(B. 


,u 


- B z v)] = 


— (B x ui 


- ByU 2 


-B zVl ). 



(10.28) 
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Since ip is an analytic function of u and v, we have, using the Cauchy-Ricmann equations, duip — d Ul ip — id U2 ip, 
and so U*dutp = U\d Ul ip — iU2 id U2 ip = Uid Ul ip + U2d U2 ip and also V*dyip = Vid Vl ip+V 2 d V2 ip. Equation (I10.13[) 
can now be written in the form 



( _ V - - A ) +eA Q + - {U x d Ul + U 2 d U2 + V x d Vl + V 2 d V2 ) 
2m \ i c J i 



(10.29) 



and after substituting (|10.16[) . we obtain for the action function a quantum Hamilton- Jacobi equation (|3.3[) 
with Hamiltonian function 



— (p 3 Aj J + eA + — — pjj + Uip Ul + U 2 pu 2 + Vip Vl + V2p V2 



(10.30) 



so that the spin coordinates have velocities u\ >2 — Ui <2 , ^1,2 = Vi,2, which agrees with Eq. (|10.21|) . Equation 
(|10.29p may be obtained from the variational principle (|10.23|) with the same Lagrangian density C as in 
Eq. (|10.24[) . but with the spin term —"fip*fSs\p there presented as (%/i)ip*(Uid Ul + U 2 d U2 + V\d Vl + V2d V2 )ip. 
The phase invariance of C leads then to the conservation law 



d f + divj + d u j u > + d U2 f* + o vl r + d V2 j V2 = 



with current 



7j?- - VW*) - — A|^| 2 , m 2 {u u U2,v 1 ,v 2 )) 

\ Aim mc J 



(10.31) 



(10.32) 



\ip\ 2 (l,v,ui,u 2 ,vi, v 2 ) 



and, therefore, to the probabilistic interpretation of |i/;| 2 and the measurement theory of section 8. 

Finally, we present the theory in "natural" coordinates Xi^^ip 011 S 3 - For that, note that u±,u 2 ,vi,v 2 in 
Eq. (110.29[) are just Cartesian coordinates in i? 4 , running from — oo to oo. To avoid confusion, rename them 
as u>i, W2, u>3, W4 and make the transformation (|10.5[) . where now ui 2 and v± j2 are real and imaginary parts of 
complex coordinates u and v, Eq. (|10.2[) . on S 3 . Let J = d(wi,W2, W3, W4)/d(x, 9, ip, r) be the jacobian matrix 
of this transformation. The spin term in the Schrodinger equation (|10.13p may then be presented, using its 
form in (|10.29j) , as 

- 7 Bs^ = - (U x d x + Ugd e + U v d v + U r d r )iP , (10.33) 

where 

(U x , U e , U v , U r ) = (Ut, U 2 ,V U V 2 ) (J^f | r=1 . (10.34) 

Since we only need this transformation on S 3 , i.e., for r = 1, Ui <2 and Vi t2 here are given by Eq. (|10.28l) where 
Mi. 2 and Ui : 2 are real and imaginary parts of u and v, Eq. (|10.2[) . Direct calculation using Mathematica then 
gives 



Uv 



U v = 



7 1 



h sin( 



(B x cos ip + By sin <p) , 



(B x sin ip — By cos ip) , 



(10.35) 



cos 8 



sine 



(-Re cos ip + By sin ip) — £> 2 



and, as expected, C/ r = 0. The Schrodinger equation now has the form 



ifi 



dip 
~dt 



1 

2 m 



eA + - (U x d x + U e d e + U v d^ 



4: 



(10.36) 
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and the Hamiltonian function in a quantum Hamilton- Jacobi equation (|3.3[) will become 

H = i ( Pj ~ c Ai ) 2 + eA ° + 2^ ftj + UxPx + Uepe + U(pPv ' (10 ' 37) 

so that x = ^ = Ue, <P = U v . 

Equation (|10.36l) may be obtained from a variational principle 5 J C dfl Ylj=i dx^dt = 0, where the measure 
dfi is given by Eq. (|10.4j) and the Lagrangian density £ by Eq. (|10.24[) with spin term —jtp*'Bstp there presented 
as (h / i)tp* (U x d x + Ugdg + U^d^) 1 ^. The measure dfl is not homogeneous — it is equal to the product of 
differentials of independent variables times the function / = sin 8. In such cases, to ensure the possibility of 
all necessary integrations by parts, the derivation of the equations of motion and conservation laws from the 
variational principle differs by using instead of the usual derivatives dk the operator djf* = (1/ f)dkf, which 
acts on any function p as di p = (1/ f)dk(fp) [Hj- In our case, when / is a function of only one variable 9, all 
derivatives except dg remain unchanged. The correct form of a conservation law that follows from the phase 
invariance of £ in a space R 3 x S 3 with an inhomogeneous integration measure d£l Ylj=i dx^ , and has the usual 
meaning and consequences there, is 

dof + divj + d x j* H 7—7; dg ( sin 9j e ) + d v f = (10.38) 
sm# 

where 

(j°J,j x ,/,^) = (h 2 , ^W*v^^vr)-^-AH 2 , \^(u x ,u e ,u v ) 

= |Vf(l,v, X ,^). 



(10.39) 



The conservation law (|10.38[) may also be obtained directly by substituting there expressions (|10.39p for 
components of the current, noticing that \ip\ 2 = e 2R , and using quantum Hamilton- Jacobi equation with 
Hamiltonian function (110.371) . 

We thus demonstrated that the analytical quantum dynamics of particles with spin can be developed by 
an extension of particle configuration space from R 3, to R 3 x S 3 . The de Broglie- Bohm- style theory of spin in 
R 3 x S 3 , in which particles are considered as a point limit of extended rigid objects, is developed in chapter 10 
of Holland's book [13]. Our spin theory uses an infinite phase space over R 3 x S 3 , and all the theory of the 
previous sections is applicable to it. In particular, particles move in R 3 x S 3 along trajectories that are well 
defined by the equations of motion, \ip\ 2 is the probability density in R 3 x S 3 , and the measurement of a spin 
component in a Stern-Gerlach experiment is a typical example of von Neumann's measurement procedure with 
discrete spectrum. 



11. Conclusion 

Let us summarize the main points of "PAQD. It is straightforward to verify that for a sum p s of a Taylor 
series (|2.6[) to satisfy PDE (|2.1[) . it is necessary and sufficient if the Taylor coefficients, or momcntums, p r a 
satisfy the ODEs (|2.5[) . where q is the velocity of an expansion point moving in the configuration space. 
Consequently, we have a simple and universal connection between evolutionary PDE (|2.ip and the dynamical 
ODEs (|2.5p . where the velocity q is still arbitrary. If, further, the Hamiltonian function H r in Eq. (|2.ip satisfies 
the Hamiltonian conditions of section 2.3, then there exists a special velocity, given by Eq. (|2.22p . which leads 
to Eq. (|2.19p . the variational principles of section 2.4, and the hierarchical Hamiltonian structure of the whole 
theory. Thus, there is a general ODE/PDE Hamiltonian formalism that may be filled with different physical 
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contents, depending on the form of a Hamiltonian function. In classical mechanics, the Hamiltonian function 
is of the first order. As was explained in section 2.7, the ODE part of the theory in this case simplifies into 
an ordinary Hamiltonian mechanics in the usual phase space. Quantum theory utilizes the second available 
option, with a Hamiltonian of a higher order and ODEs residing in an infinite phase space. More specifically, 
it appears that in a nonrelativistic domain for spinless particles, nature builds quantum theories by the recipe 
of section 4, so that any theory of this kind is defined by Eq. (|4.3p with some Lagrangian function, quadratic 
in velocity and such that the corresponding Hamiltonian function satisfies Hamiltonian condition (|2.24l) . The 
theory will then automatically exhibit the superposition principle, path-integral representation, wave-particle 
duality (which is shown to be possible only in the infinite phase space), and the classical limit, described by 
Hamilton- Jacobi PDE and Hamilton ODEs with a Hamiltonian, corresponding to the Lagrangian function 
in Eq. (|4.3[) . Also, the resulting Schrodinger equation will be automatically obtainable from a variational 
principle, so its symmetries will lead to corresponding conservation laws. Since Hamiltonian functions in our 
theory depend only on derivatives of unknown functions, they are automatically invariant with respect to 
shifts of these functions by arbitrary constants. This symmetry leads to a current conservation, a current 
being defined with the correct velocity (|2.22j) . and to invariance of a measure |-0| 2 dy with respect to equations 
of motion. This invariance leads then to the probability density p = \ip\ 2 in the same way as invariance of the 
Liouville measure leads to the microcanonical distribution in classical statistics, the difference in probability 
densities resulting from different forms of equations of motion. The probabilistic interpretation of the wave 
function is, therefore, deduced in "PAQD rather than being postulated. The multiparticle generalization of the 
theory leads to the standard picture of quantum particles in a classical macroscopic environment, and being 
applied to specially constructed apparatuses, to the quantum theory of measurement. The measurements of 
classical quantities that may be used as parts of particle-apparatus interaction Hamiltonians, appear to have the 
desired features in this theory. On the other hand, the measurement and observation of nonlocal momcntums, 
such as PxiVj , where i ^ j are particle indices, is impossible, because Hamiltonians built by the rules of section 
4 cannot contain such terms. The presence of such momentums makes the whole multiparticle theory nonlocal, 
and explains the mechanism of nonlocal correlations. On the other hand, their nonobservability prevents using 
them for the transmission of superluminal signals. As was discussed in section 9, in spite of being nonlocal, the 
relativistic version of the theory seems to be presentable in a Lorentz invariant and even generally covariant 
way. Finally, the theory of particles with nonzero spin resides in configuration space that includes, along with 
particle's space coordinates, its internal SU(2) degrees of freedom. 

In Table 1, different aspects of "PAQD and QM are compared in a self-explanatory form. The comparison 
clearly demonstrates that in spite of experimental agreement, the two theories are evidently different and draw 
different pictures of the physical reality. Several additional remarks may be useful. First, as was discussed in 
sections 2.5 and 3, the actual solution of the equations of motion may be obtained by the generalized Jacobi 
method from a known wave function. In the spinless case, the particles will then move along Bohmian trajec- 
tories. The equations of motion, however, determine the motion completely and unambiguously themselves, 
and so a technique of their direct solution, without any use of a wave function, should be possible to develop if 
desirable. Second, the theory is formulated in terms of momentums, and as was just discussed not all of them 
are observable. We note, however, that all coordinates and momentums that are observable (for example in 
such devices as bubble chambers) are reflected in the theory, and nonobservable momentums are nonobservable 
not because they are postulated to be such, but because this is a property of an observation/measurement 
procedure that follows from the basic equations of the theory. The situation here should be compared with 
the one in QM, which is formulated entirely in terms of nonobservable wave functions, and brings in the ob- 
servable quantities (coordinates and momentums) only through the interpretational part of the theory. Third, 
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PAQD 


QM 


Relation of mathematical theory to experiment 


Straightforward 


Needs interpretation 


Laws of nature 


Deterministic 


Indeterministic 


Particles behavior 


Always particle-like 


Complementary 


Wave-particle duality 


Matncmatically derived 


Verbally described 


INon-statistical interpretation ot wave lunction 


Exists 


Does not exist 


Statistical interpretation of wave function 


Mathematically derived 


Postulated 


Measurement problem 


Does not exist 


Unsolved 


Nonlocality vs. Lorentz invariance conflict 


Cleared 


Unresolved 


Classical limit 


Direct 


Indirect 


Geometric picture of classical mechanics 


Generalized 


Lost 



Table 1: The comparison of "PAQD and QM. 



the theory of von Neumann's measurements, presented in section 8, exhibits all the properties attributed to 
the measurement procedure in standard QM, in particular such measurements must satisfy the uncertainty 
relations. However, contrary to QM, "PAQD gives detailed description of the behavior of both observed system 
and observing apparatus before, during, and after their interaction. In this situation, it is not unthinkable to 
speculate that new "measurement-like" procedures may be found, which will provide more information than is 
permitted by the uncertainty principle, or will generate experimental situations (especially when mesoscopic 
objects are involved) in which QM (but not PAQD) fails to give unambiguous predictions. The investigation 
of such possibilities lies, however, outside of the scope of the present work. 

To summarize, the mathematical theory developed in section 2 allows to give a simple description of 
quantum phenomena as resulting from generalized Hamiltonian motion of particles. The present theory does 
not suffer from the shortcomings discussed in the Introduction. It appears especially important, that although 
the specific form of the theory is completely determined by the Schrodinger equation, it nevertheless allows 
to simultaneously deduce the statistical interpretation, which in existing quantum theory is described by the 
separate and independent axioms. We conclude, that it seems not unreasonable to believe that PAQD may 
indeed provide the true and fundamental description of nature. 
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Appendix: Statistical distribution and entropy growth in classical 
statistics 

To facilitate the comparison in section 6 of statistical distribution in PAQD and in the classical theory, here 
we review the basics of classical equilibrium and nonequilibrium statistics in a convenient for this comparison 



59 



form. A classical system is represented by a point in a phase space P of dimension 2n p n s , where n p is the 
number of particles and n s is the dimension of the physical space. Denote a general point of P by 7. Invariant 
measure in P is a Liouville measure Y[7=i d ns rid ns pi, and we denote it by dT. Besides the description in 
a phase space P, or the space of microstates that we will also call "full description space" (FDS), classical 
statistics uses a much cruder description of systems by sets of kinetic or thermodynamic parameters forming 
"reduced description space" (RDS) or the space of macrostates of dimension ur <C 2n p n s . The coordinates 
Ri of this space are the values of some functions ^(7), i = 1, . . . ,ur of the microstate, and the macrostate 
is considered fully specified by the known values of all the Ri, i.e., by the set R = {Ri,i = 1, . . . ,ur}. In 
other words, phase space P is broken into subspaces Pr that correspond to small cells in RDS: Pr = {7 € 
P : Ri < 7*1(7) — Ri + * = 1 , - - - , tir}, and a crude description of a system at any time t is given by the 
corresponding R(t), i.e., by specifying a subspace that the point 7 resides in at this moment. Different points 
of Pr represent then different microstates, compatible with the same macrostate R, so that every time the 
system is prepared in this macrostate, its microstate will be represented by some random point 7 G Pr. As 
it was with a quantum particle, in the mathematical limit of an infinite number of such preparations, these 
points will form an ensemble Ar with probability density pr that maximizes the corresponding Gibbs entropy 
S G R \p). Repeating the steps that led from (I6.1[) to (|6.3p . but using the Liouville measure this time, it is easy 
to show that 

S ( g\p) = - I p\n P dT. (Al) 

JPr 

We call the density pr that maximizes this Sq (p) a microcanonical density, and corresponding ensemble Ar 
a microcanonical ensemble. Let us show that the microcanonical density is constant: 

Pr = const = — — , (A2) 

where T(Pr) = J p ^ dT is a phase volume of Pr. Indeed, from pr — const and In a; < x — 1, we have for any 
other normalized probability density p in Pr: 

Pr\i\ PRdT = — I pin p dT p In — dT 

Pr ' JPr JPr P 

> -J plnpdT + J p(l-^dT (A3) 
= — / plnpdT . 

JPr 

Due to the constancy of pr, the averaging over Ar, i.e., over a set of systems randomly created in a 
macrostate R, reduces to the averaging over Pr with the measure dr. Consequently, the need for introducing 
ensembles Ar disappears; instead, we will use as ensembles the corresponding subsets Pr of FDS. 

The maximum possible value of the Gibbs entropy on Pr will be called the Boltzmann entropy of R and 
denoted as Sb{R)- We have, obviously, 

S B (R) = max S ( G R) (p) = S { G R) (p R ) = -lnp R = lnT(P R ) . (A4) 
{p} 

While Gibbs entropy is a function of the probability density p in FDS, Boltzmann entropy is a function of a set 
R that belongs to a RDS and describes a macroscopic state of a system. The corresponding functions rj may 
be, for example, the particle, energy, and momentum densities in small cells, covering the volume of a system 
or, for another example, the densities of particles in the cells that cover a six-dimensional one-particle phase 
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space (i.e., essentially, the values of the Boltzmaim distribution function in different points of this space). 
The fact that the Boltzmann entropy Sb depends on the macroscopic state R of a system, rather than on its 
microscopic state 7, is a manifestation of its "anthropomorphic" nature, as was discussed by Jaynes by 
its very thermodynamic definition, the difference of a system's entropy between different states depends on 
which parameters were held fixed and which were allowed to change during the transition from one state to 
the other. 

The relations (|A2[) and (|A4I) are based only on the N — > 00 limit (see Eq. (|6.4I) ) . They are, therefore, 
universally applicable to any system and any functions r i: even if they don't have a macroscopic charac- 
ter. In macroscopic systems, however, functions r.i may be selected in a way that allows them to possess 
additional important properties |30] . Namely, functions give a crude description of a system; as was just 
discussed, usually they represent the properties of particles in small cells in coordinate or phase spaces. Then 
in macroscopic systems, with their enormous quantity of particles, the size of these cells may be chosen in 
such a way that on one hand they are sufficiently small, so that the parts of the system inside them appear 
homogeneous and further division does not lead to more refined description, while on the other hand they are 
large enough to still contain a macroscopic number of particles. The behavior of all physical values that are 
additive with respect to contributions of separate particles or small groups of particles, will then be regulated 
by a central limit theorem: in a typical microstate 7 S Pr, these values, as well as their time derivatives, 
will be close to the corresponding averages over Pr with negligible dispersion. The averaging over Pr serves, 
therefore, as a method of calculation of typical, i.e., observable, magnitudes of thermodynamic values, with 
the microcanonical ensemble often being replaced by a canonical one for calculational convenience. 

The existence of such intermediate scale of description refinement, and of corresponding typical behavior, 
practically identical to average one, is the first important property of macroscopic systems. Another property 
is related to the character of their time evolution. Namely, let at time to the system be described by parameter 
set R(to), so that its microstate belongs to an ensemble (i.e., subspace of FDS) PR(t )- We let all states of this 
ensemble evolve until the time t\ > to and denote the resulting ensemble PRU \(ti). As was just discussed, 
typical values of the functions r"j in this ensemble will be close, with negligible dispersion, to corresponding 
averages over it, which we will denote by Ri(ti). The typical, i.e., observable, macrostate R(t\) for t\ > to is 
obtained, therefore, by direct averaging of the equations of motion over Pr(j ). We refer to the corresponding- 
averaged equations as generalized kinetic equations. Let us now compare the ensemble PR(t )(ti) with a 
microcanonical ensemble Pr(h) that corresponds to the values of the parameters R, observed at time t\. 
Neglecting extremely rare nontypical microstates, we can say that the values of R in all states 7 € PR(t )(ti) 
are equal to R(t\), so that PR(t ){ti) is a subset of Pr^), PR(t ){ti) C PR(t t )- We need to elaborate in what 
sense the microstates of PRtt \{tx) that do not belong to Pr^) are "extremely rare." Note that phase volumes 
that are essential for our conclusions are such that their logarithms are extensive, i.e., scale proportionally 
to the number of particles n p in a system. Consequently, we are only interested in the logarithms of these 
volumes per particle in the limit as n p —> 00. Then the inclusion PR(t Q ){ti) C Pr^) should be understood as 

lim ^fl n r(P R(to) (t 1 ))-lnr(P R{to) (t 1 )r\P R( t 1 )))=0, (A5) 

Yl-p rOO Tip V / 

so that the part of Pft( to )(ii) that falls outside of P R (ti) is inessential in this sense of logarithm per particle. 
This consideration also explains why we should not care about the exact values of the Si in the definition of 
Pr — they are inessential in the same sense. 

We are especially interested in the situation when at time to the system was in a nonequilibrium state. The 
states of the ensemble P.R(t )(ti) will not be typical for Pru^ m this case, for they will have nontypical for 
Pfj( tl ) correlations. These correlations will manifest themselves under time inversion: the states of Pfj( to )(fi) 
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will return back to the macrostate R{to), which is further from equilibrium than R(ti), while the typical states 
of Pfl(ti) will approach equilibrium (modulus tiny thermodynamic fluctuations) with deviation of time from t\ 
in both directions. Also, if we wait for the whole Poincare cycle to pass, then we will observe another manifes- 
tation: the states of P R ^(ti) will deviate off equilibrium as far as to R(to), while the typical states of PrUx) 
will deviate less — up to R(t\). It happens, however, and this is the second important property of macro- 
scopic systems, that the correlations, which are different in P R ^ to -\(ti) and Pfj( tl ), are n-particle correlations 
with macroscopically large n ("large-n correlations"), while "small-n correlations" (i.e., one-particle densities 
and correlations between a small number of particles) in PR(t )(ti) and Pr(h) are practically the same. At the 
end of the Poincare cycle nontypical large-n correlations will conspire to coherently affect small-n ones and 
drive the system into an abnormally (for Pr^^) nonequilibrium state R(to), but before that the influence of 
large-n correlations on small-n ones may be neglected. But it is only these small-n correlations, and not large-n 
ones, that contribute to observable and measurable physical values of interest. Consequently, in the normal 
physical experiment, when time goes only forward, but not as far forward as for the length of the Poincare 
cycle, nontypical for PR(ti) large-n correlations in PR(t ){t\) do not manifest themselves, and we can regard the 
states of PR(t )(ti) as typical states of Pr^)- In other words, the origin of the current microstate 7(^1) of the 
system, reflected in the ensemble -PR(t„)(ti), becomes irrelevant, the only important question being in which 
subspace Pr of P point 7 resides now. This means that for every time t > to, ensemble PR( to )(t), representing 
the system, can (and must, for that makes the calculations much simpler) be replaced by a microcanonical 
ensemble Pmt), where R(t) is obtained from initial R(to) by solving the generalized kinetic equations. Now 
from the inclusion PRU \(t) C PR(t) an d the invariance of the Liouville measure, we have for t > to 

r{P a( t )) = r(p R(to) (t)) < r(p m ) (A6) 

and so <Sb(-R(<)) > SB(R{to)), i.e., Boltzmann entropy never decreases and achieves its maximum in equilib- 
rium, when the system's macrostate R ceases to change. 
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